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Abstract 

The  primary  goal  of  this  research  was  to  determine  the  feasibility  of  air-to-air 
optical  communication  systems  using  coherent  detection  in  a  turbulent  environment. 
Secondary  goals  were  to  determine  (1)  the  signal- to-noise  ratio  (SNR)  degradation 

f 

due  to  turbulence  and  (2)  the  effect  of  varying  the  wavelength  (1.5,  830,  904  nm), 
altitude  (5000-12500  m),  path  length  (40-160  km)  and  aperture  side  length  (1-15  cm) 
on  the  turbulence  affected  SNR. 

Cl  The  research  was  conducted  under  the  following  assumptions:  (1)  ideal  trans¬ 
mitter  and  receiver,  (2)  equally  sized,  uniformly  weighted,  square  transmitter  and 
receiver  apertures,  (3)  ideal  tracking  by  the  receiver,  and  (4)  zero  losses  due  to  beam 
steering,  beam  spreading  and  scintillation/^) 

— It  was  shown  that  the  SNR  efficiency  (defined  as  the  ratio  of  the  turbulence 
affected  SNR  to  the  non-turbulence  affected  SNR)  ranged  from  0.16  to  67.2  percent. 
The  67.2  percent  efficiency  was  achieved  for  a  wavelength  of  904  nm,  path  length  of 
40  km,  altitude  of  5000  m  and  aperture  side  length  of  1  cm.  f _ .... 


The  feasibility  analysis  showed  that,  depending  on  the  transmitter  power  and 
data  rate,  path  lengths  in  excess  of  160  km  could  be  achieved  for  bit  error  rates  of 
both  10-5  and  10-10. 
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CHARACTERIZATION  OF  AN 
AIR-TO-AIR  OPTICAL  HETERODYNE 
COMMUNICATION  SYSTEM 


I.  Introduction 

1,1  Background 

In  air-to-air  applications,  optical  communication  systems  have  shown  promise 
as  a  means  of  covert  communications  with  high  data  rate  capability.  Experiments, 
such  as  HAVE  LACE  (Laser  Airborne  Communications  Experiment)  conducted  by 
the  United  States  Air  Force,  have  shown  that  air-to-air  optical  communication  sys¬ 
tems  employing  incoherent  (direct)  detection  are  feasible  with  acceptable  signal-to- 
noise  ratio  (SNR),  and  thus  probability  of  error,  performance.  However,  the  SNR 
achievable  in  incoherent  detection  systems  is  quantum  limited;  that  is,  there  is  a 
“...  minimum  pulse  energy  Emin  required  to  maintain  a  given  bit  error  rate  (BER) 
which  any  practical  receiver  must  satisfy  ...”  (33:356). 

Coherent  (heterodyne)  detection  systems,  on  the  other  hand,  offer  a  significant 
improvement  in  SNR  performance  on  the  order  of  10-20  decibels  (dB)  (25,  26,  32). 
Greater  SNR  performance,  in  turn,  leads  to  higher  data  rate  capability  and  better 
probability  of  error  (or  BER)  performance. 

Both  incoherent  and  coherent  detection  systems,  however,  are  adversely  af¬ 
fected  by  the  effects  of  atmospheric  turbulence  on  the  propagating  optical  wave. 
The  resulting  amplitude  fluctuations  of  the  optical  wave  due  to  turbulence  affect 
both  incoherent  and  coherent  detection.  The  effect  is  reduced  power  (or  SNR)  at 
the  receiving  aperture.  In  addition,  coherent  detection  is  further  affected  by  the 
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optical  wave’s  phase  fluctuations  due  to  turbulence.  In  coherent  detection,  a  locally 
generated  wave  (local  oscillator)  is  combined  with  the  received  wave.  The  portions 
of  the  received  wave  that  are  out  of  phase,  due  to  phase  fluctuations,  with  the  detec¬ 
tor’s  local  oscillator  will  combine  destructively  resulting  in  an  even  further  reduction 
in  the  received  SNR. 

The  Communications  Technology  Branch  of  the  Wright  Research  and  Develop¬ 
ment  Center  (WRDC),  Wright-Patterson  AFB,  OH  is  presently  sponsoring  research 
to  determine  the  effects  of  atmospheric  turbulence  on  optical  communication  systems 
employing  coherent  detection  systems. 

1.2  Problem 

The  basic  problem  facing  WRDC  is  to  determine  the  feasibility  of  air-to-air 
optical  communication  systems  using  coherent  detection  in  a  turbulent  environment. 
This  study  will  characterize  the  performance  of  an  optical  heterodyne  communication 
link  between  aircraft  in  the  presence  of  atmospheric  turbulence.  If  optical  heterodyne 
communication  systems  prove  to  be  feasible  in  air-to-air  applications,  these  systems 
will  provide  a  means  of  covert  communications  with  higher  data  rate  capability 
than  conventional  radio  frequency  and  microwave  communication  systems  currently 
employed  in  aircraft,  as  well  as  incoherent  optical  communication  systems. 

1.3  Definitions 

Before  the  scope  and  assumptions  of  this  study  can  be  presented,  several  terms 
and  concepts  must  be  defined. 

Power  spectral  density  of  turbulence  (  $„(«;)):  The  power  spectral  density 
(PSD)  of  turbulence  is  a  measure  of  the  number  of  eddies  of  different  sizes 
making  up  the  turbulence  (18:388).  The  PSD  is  designated  as  $„(«)  where  k 
is  the  wave  number  (k  =  27r/L  where  L  is  the  eddy  size).  Eddies  are  packets 
of  air  in  which  the  refractive  index  of  each  of  the  packets  varies  with  respect 
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to  the  average  index  of  refraction  (10,  18).  The  spectrum,  $n(/c),  consists  of 
three  regions  determined  by  the  inner  and  outer  scales  of  the  turbulence  which 
are  directly  related  to  the  eddy  sizes.  In  the  region  consisting  of  eddy  sizes 
greater  than  the  outer  scale,  $n(/c)  may  be  described  by  either  the  von  Karman 
spectrum  or  the  Kolmogorov  spectrum.  In  the  region  consisting  of  eddy  sizes 
between  the  inner  and  outer  scales,  $n(/c)  is  described  by  the  inertial  subrange 
of  the  Kolmogorov  spectrum.  Finally,  in  the  region  consisting  of  eddy  sizes  less 
than  the  inner  scale,  $„(ac)  is  again  described  by  the  Kolmogorov  spectrum. 
Chapter  II  contains  a  complete  discussion  of  $n(/c)  as  well  as  an  illustration  of 
the  regions  (Figure  2)  (18:388-389). 

Receiver  fundamental  limits:  Receiving  optical  systems  axe  fundamentally  lim¬ 
ited  in  performance  capability  due  to  the  basic  laws  of  nature  (that  is,  quantum 
mechanics)  (30:262-263). 

Isotropic  and  homogeneous  turbulence:  When  the  refractive  index  fluctuations 
of  the  turbulence  have  a  spherically  symmetric  autocorrelation  function,  the 
turbulence  is  statistically  isotropic.  Similarly,  when  the  refractive  index  fluc¬ 
tuations  are  stationary  from  a  three-dimensional  spatial  standpoint,  the  tur¬ 
bulence  is  statistically  homogeneous  (18:386). 

Frozen  atmosphere:  Under  the  assumption  of  a  frozen  atmosphere  (or  frozen  tur¬ 
bulence),  the  atmospheric  refractive  index  fluctuations  do  not  vary  with  time, 
but  rather  drift  across  the  channel  with  constant  velocity  due  to  wind  (18:386). 

Refractive  index  fluctuations  structure  constant  (C2):  The  refractive  index 
fluctuations  structure  constant,  C2  (m~2/3),  is  a  measure  of  the  strength  of  the 
fluctuations  and  depends  on  both  the  atmospheric  conditions  and  the  altitude 
of  the  turbulence. 

Mutual  Coherence  Function  (MCF):  TheMCF  describes  the  loss  of  coherence 
of  a  coherently  transmitted  wave  (29:482).  Of  special  interest  to  this  research 
effort  is  the  loss  of  coherence  caused  by  propagation  through  turbulence. 
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Cross-Correlation  Function  (if(p,r,z)):  H(p,r,z)  is  the  two  point  correlation 
function.  For  finite  optical  beam  propagation,  all  of  the  propagation  properties 
needed  to  determine  the  MCF  are  contained  in  J/(£,  r,  z)  (28:2154).  The  terms 
g  and  r  represent  points  in  the  transmitter  and  receiver  planes,  respectively; 
while  z  represents  the  distance  between  the  two  planes. 

1.4  Scope 

This  study  investigates  the  feasibility  of  an  air-to-air  optical  heterodyne  com¬ 
munication  link  from  a  probability  of  error,  or  BER,  performance  standpoint.  The 
communication  system  parameters  to  be  varied  are  transmitter/receiver  aperture 
sizes,  aircraft  separation,  aircraft  altitudes  and  data  rate. 

The  turbulence  is  modeled  based  on  Gaussian  cross-correlation  statistics  for 
the  turbulence  induced  phase  fluctuations.  The  turbulence  strength,  C£,  described 
in  Chapter  II,  is  modeled  using  the  Hufnagel  model  (23:6-14). 

This  study  only  considers  the  turbulence  losses  due  to  phase  fluctuations  re¬ 
sulting  from  propagation  through  turbulence.  Losses  due  to  beam  spreading,  beam 
steering  and  scintillation  axe  not  addressed.  This  study  also  does  not  consider  aircraft 
skin  turbulence  caused  by  aircraft  motion. 

1.5  Standards 

For  this  research  effort,  a  BER  of  10“5  is  considered  acceptable  for  feasibil¬ 
ity.  This  value  is  commensurate  with  similar  research  and  experimental  efforts.  For 
example,  HAVE  LACE  required  a  BER  of  10“6  at  160  km  (22:6),  while  Kanavos 
required  a  BER  of  10"5  during  a  study  of  air-to-air  direct  detection  optical  com- 
munciation  systems  (24:58,63). 

1.6  Assumptions 

This  study  is  conducted  with  the  following  assumptions: 
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1.  An  ideal  transmitter  and  receiver. 


2.  Equal  transmitter  and  receiver  dimensions. 

3.  Ideal  tracking  by  the  receiver. 

4.  Locally  isotropic  and  locally  homogeneous  turbulence.  Note  that,  ac¬ 
cording  to  Davis,  turbulence  consisting  of  eddy  sizes  between  the  inner 
and  outer  scales  of  turbulence  is  both  isotropic  and  homogeneous  (9:140). 
However,  as  discussed  by  Goodman,  turbulence  consisting  of  eddy  sizes 
greater  than  the  outer  scale  is  both  nonisotropic  and  nonhomogeneous, 
but  the  refractive  index  fluctuations  may  still  possess  an  isotropic  and 
homogeneous  structure  function  (18:391). 

5.  “Frozen”  atmospheric  turbulence.  This  assumption,  also  known  as  Tay¬ 
lor’s  hypothesis,  follows  from  the  fact  that  the  “fluctuation  time”  of  the 
turbulence  is  much  longer  than  the  time  required  for  the  optical  wave  to 
propagate  through  the  atmospheric  communication  channel  (18:386). 

6.  A  Gaussian  probability  density  function  (PDF)  is  used  to  characterize 
the  statistics  of  the  phase  fluctuations  resulting  from  the  optical  wave 
propagating  through  the  atmospheric  turbulence.  This  characterization 
is  discussed  further  in  Chapter  II. 

7.  Square  transmitter  and  receiver  apertures.  The  SNR  equation,  developed 
in  Chapter  III,  can  be  reduced  to  a  four-fold  integral  using  circular  aper¬ 
tures,  and  2  two-fold  integrals  using  square  apertures.  Due  to  the  limited 
processing  capability  of  even  today’s  computers,  current  numerical  inte¬ 
gration  methods  perform  best  on  only  three-fold  or  less  integrals,  thus 
discouraging  their  use  for  the  circular  aperture  SNR  equation. 

8.  A  Gaussian  cross-correlation  function,  H (p,  r,  z).  The  use  of  the  Gaussian 
cross-correlation  function  allows  the  SNR  equation  to  be  reduced  to  the 
2  two-fold  integrals  for  the  square  apertures  discussed  above. 
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9.  Zero  losses  due  to  beam  steering,  beam  spreading  and  scintillation.  This 
was  discussed  in  Section  1.4. 

1.7  General  Approach 

To  determine  the  feasibility  of  an  optical  heterodyne  communication  link  in 
the  presence  of  atmospheric  turbulence,  the  following  approach  is  taken: 

1.  Investigate  the  relationship  between  aperture  dimensions  and  SNR. 

2.  Develop  a  theoretical  model  of  the  SNR  both  with  and  without  the  effects 
of  atmospheric  turbulence  considered. 

3.  Determine  a  SNR  efficiency  factor,  a,  which  is  the  ratio  of  the  turbulence 
affected  SNR  to  the  non-turbulence  affected  SNR. 

4.  Once  the  SNR  model  is  developed,  calculate  the  SNRs  (with  and  without 
turbulence  effects  considered)  while  varying  the  following  parameters: 

(a)  Transmitting  and  receiving  aperture  sizes.  The  side  lengths  (of 
the  square  apertures)  considered  range  between  1.0  and  15.0  cm. 

(b)  Aircraft  separation.  Separations,  or  path  lengths,  considered 
range  from  40  to  160  km. 

(c)  Aircraft  altitude.  Altitudes  considered-range  from  5000  to  12500  m. 

(d)  Wavelength.  Wavelengths  considered  are  1.5,  830,  and  904  nm. 

(e)  Aircraft  velocity.  Velocities  considered  range  from  141.0  to  244.6 
m/s,  although  the  majority  of  the  analysis  uses  218.0  m/s. 

(f)  Data  rate.  Data  rates  considered  are  20000,  40000,  100000  and 
1000000  bps. 

5.  Calculate  the  BERs  from  the  associated  SNRs. 

6.  Evaluate  the  BERs  to  determine  the  feasibility  of  optical  heterodyne  com¬ 
munication  systems  in  the  presence  of  atmospheric  turbulence. 
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1.8  Sequence  of  Presentation 

This  research  report  is  organized  as  follows: 

1.  Chapter  II,  the  Literature  Review,  presents  a  summary  of  the  current 
knowledge  in  the  areas  of  optical  communication  systems  and  coherent 
detection  systems  in  the  presence  of  atmospheric  turbulence. 

2.  Chapter  III  presents  the  theory  and  methodology  applied  in  developing 
the  models  outlined  in  the  General  Approach. 

3.  Chapter  IV  contains  the  analysis  and  evaluation  of  the  calculated  SNRs 
and  BERs.  From  the  BERs,  feasibility  is  determined. 

4.  Chapter  V  contains  the  conclusions  of  the  feasibility  evaluation  and  rec¬ 
ommendations  for  further  research. 

In  addition,  this  report  contains  several  appendices. 

1.  Appendix  A  contains  a  discussion  of  the  different  nomenclature  and  math¬ 
ematical  expressions  used  to  describe  the  phase  coherence  length. 

2.  Appendix  B  contains  the  development  of  the  SNR  for  circular  apertures. 

3.  Appendix  C  contains  the  development  of  the  computer  code  which  models 
the  SNRs.  This  appendix  also  contains  a  copy  of  the  computer  code. 
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II.  Literature  Review 


2.1  Introduction 

Since  the  development  of  the  laser,  there  has  been  a  resurgence  of  interest  in 
the  use  of  optical  communication  systems.  Guided  optical  (or  fiber  optic)  systems, 
for  example,  are  being  used  extensively  in  telephone  and  data  communications.  In 
the  area  of  unguided  systems,  eaxth-to-space  and  space-to-space  optical  communi¬ 
cation  systems  for  high  data  rate  transfer  have  been  developed  (16:22-3).  Also, 
air-to-ground  (Air  Force  405B  project)  and  air-to-air  (HAVE  LACE)  optical  com¬ 
munication  experiments  have  been  conducted  by  the  United  States  Air  Force  (10:1). 
This  research  effort  focuses  on  unguided  optical  systems  operating  in  the  atmosphere. 

Both  the  Air  Force  HAVE  LACE  program  and  theoretical  research  have  shown 
that  optical  communications  can  be  conducted  through  the  earth’s  atmosphere  with 
acceptable  signal-to-noise  ratio  (SNR)  performance  (9,  10,  11,  16,  24).  The  use  of 
optical  communications  systems  in  the  role  of  air-to-air  communications  between 
aircraft  offers  many  advantages,  as  well  as  several  disadvantages,  with  respect  to 
radio  frequency  (RF)  and  microwave  communications.  Table  1  lists  several  of  these 
advantages  and  disadvantages. 

Three  of  the  advantages  listed  are  of  particular  interest  to  the  Air  Force: 

1.  Small  receiver/transmitter  apertures:  small  receiver  and  transmitter  aper¬ 
ture  sizes  are  desirable  due  to  the  weight  and  size  limitations  aboard 
strategic  and  tactical  aircraft. 

2.  High  data  rate  capability:  Feldmann  states  that  communications  systems 
performing  “Strategic  Data  Exchange”  (for  example,  intelligence  collec¬ 
tion  or  command  and  control  missions)  will  be  required  to  transmit  data 
at  rates  ranging  from  19.2  kilobits  per  second  up  to  several  megabits  per 
second  (10:1-2). 
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Table  1.  Advantages  and  Disadvantages  of  Optical  Communications  Over  Radio 
Frequency  and  Microwave  Communications  in  the  Atmosphere  (10, 19, 24) 


Advantages 

Disadvantages 

-  Higher  Received  Power  to 

Transmitted  Power  Ratio 

-  Small  Receiver/Transmitter 

Apertures 

-  High  Data  Rate  Capability 

-  Low  Probability  of  Intercept/ 

Increased  Jam  Resistance 

-  Better  Directional  Resolution 

With  Smaller  Apertures 

-  Exploitation  of  the  Unused 

Portion  of  the 
Electromagnetic  Spectrum 

-  No  Communication  License 

Required 

-  Atmospheric  Effects 

-  Less  Suited  to  Broadcasting 

Due  to  Narrow  Beam 

-  Accurate  Pointing  and 

Tracking  Required 

-  Higher  Level  of  Noise  in 

Receiver  Due  to  Quantum 
Limited  Nature  of  Detection 
Process 

-Low  Generation  Efficiency 
in  Optical  Carrier 

-  Possible  Hazard  (Due  to 

Laser  Beam) 

3.  Low  probability  of  intercept/  increased  jam  resistance:  in  both  strategic 
and  tactical  environments,  covert  communications  may  be  necessary  to 
accomplish  the  mission,  thus  the  need  for  communications  systems  with 
low  probability  of  intercept.  In  the  command  and  control  role,  communi¬ 
cation  is  paramount,  thus  calling  for  jam  resistant  systems. 

On  the  other  hand,  the  primary  disadvantage  listed  in  Table  1,  with  respect 
to  air-to-air  communications  between  aircraft,  is  the  atmospheric  effects.  The  atmo¬ 
spheric  elfects  on  an  optical  communication  link  may  severely  limit  the  range  and 
the  probability  of  error  performance  of  an  optical  communication  system. 

The  following  sections  of  this  literature  review  address  the  effects  of  the  atmo¬ 
sphere  on  an  optical  communication  system.  Section  2.2  covers  atmospheric  turbu¬ 
lence  to  include  observable  effects  on  an  optical  wave  propagating  through  turbu¬ 
lence.  Section  2.3  addresses  detection  of  an  optical  wave,  and  the  effects  of  turbulence 
on  the  detection  process.  Section  2.4  then  discusses  current  research  efforts  in  model- 
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ing  both  the  MCF  and  SNR  with  respect  to  coherent  detection.  Section  2.5  addresses 
BER  with  respect  to  direct  detection  systems.  Finally,  Section  2.6  presents  a  short 
summary. 

2.2  Atmospheric  Turbulence 

Atmospheric  turbulence  can  be  viewed  as  packets  of  air,  called  eddies,  in  which 
the  refractive  index  of  each  of  the  eddies  varies  with  respect  to  the  average  index 
of  refraction  (10,  18).  Figure  1  illustrates  an  optical  beam  propagating  through 
the  eddies.  As  the  optical  beam  propagates,  it  experiences  both  phase  and  ampli¬ 
tude  fluctuations.  The  phase  and  amplitude  fluctuations  result  in  several  observable 
effects  on  the  propagating  wave  front  (9:139): 

1 .  Beam  Steering:  the  entire  beam  may  be  deviated  from  the  line  of  sight, 
giving  rise  to  a  loss  of  power  at  the  receiving  aperture. 

2.  Beam  Spreading:  small-angle  scattering  by  the  inhomogeneities  spreads 
the  signed  energy  over  a  large  region,  resulting  in  a  reduced  energy  density 
at  the  receiver. 

3.  Scintillation:  destructive  interference  within  the  beam  produces  local 
fluctuations  in  the  amplitude,  so  that  over  the  beam  there  appear  areas 
brighter  and  darker  when  compared  to  the  average.  This  leads  to  fluctu¬ 
ations  in  received  power  which  are  strongly  dependent  on  the  dimensions 
of  the  receiver  aperture.  The  time  variations  of  these  fluctuations  must 
be  considered  in  amplitude  modulation  (AM)  signals. 

4.  Phase  Fluctuations:  the  time  variation  in  phase  produces  a  spurious  mod¬ 
ulation  which  interferes  with  the  signal  modulation,  and  the  spatial  vari¬ 
ation  tends  to  destroy  the  signal  modulation  in  AM  signals. 
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Tx  -  Transmitter 
Rx  -  Receiver 

-  Without  Turbulence 

-  With  Turbulence 
^  -  eddies 


With  respect  to  coherent  detection  systems,  beam  steering  and  phase  fluctua¬ 
tions  are  of  primary  interest.  This  research  effort  concentrates  on  the  latter,  phase 
fluctuations,  caused  by  atmospheric  turbulence. 

2.2.1  Atmospheric  Turbulence  Description  The  refractive  index  of  air 
at  optical  frequencies  is  dependent  on  the  optical  beam  wavelength,  atmospheric 
pressure,  and  temperature  (8,  18).  Many  authors  define  n\  as  a  measure  of  the 
deviation  of  the  refractive  index  from  the  free  space  value.  The  free  space  value  is 
approximately  equal  to  unity  (8, 18). 

According  to  Goodman,  the  variations  on  the  refractive  index  due  to  pressure 
are  small  and  can  be  ignored  (18:388).  Therefore,  the  refractive  index  variations  are 
primarily  due  to  temperature  variations.  How  the  temperature  variations  result  in 
refractive  index  inhomogeneities  is  best  explained  by  Goodman: 


The  random  fluctuations  n\  of  the  refractive  index  are  caused  predomi¬ 
nantly  by  random  microstructure  in  the  spatial  distribution  of  temper¬ 
ature.  The  origin  of  this  microstructure  lies  in  extremely  large  scale 
temperature  inhomogeneities  caused  by  differential  heating  of  different 
portions  of  the  Earth’s  surface  by  the  sun.  These  large-scale  tempera¬ 
ture  inhomogeneities,  in  turn,  cause  large-scale  refractive  index  inhomo¬ 
geneities,  which  are  eventually  broken  up  by  turbulent  wind  flow  and 
convection,  spreading  the  scale  of  the  inhomogeneities  to  smaller  and 
smaller  sizes.  (18:388) 


The  refractive  index  inhomogeneities  affect  an  optical  wave  in  two  distinct 
ways:  distortion  of  the  phase  front  such  that  the  isophase  surface  of  the  wave  front 
is  no  longer  plane,  and  fluctuations  in  the  amplitude  or  intensity  across  the  wave 
front  (3,  9, 11, 12).  Recall,  it  is  the  phase  and  amplitude  fluctuations  that  cause  the 
observable  effects  listed  in  the  previous  section. 

2.2.2  Atmospheric  Turbulence  Models  Because  of  the  random  nature  of 
the  refractive  index  fluctuations,  a  deterministic  model  of  the  atmospheric  effects  is 
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not  possible.  That  is,  the  refractive  index  eddies  vary  randomly  not  only  in  size  and 
location,  as  illustrated  in  Figure  1,  but  also  with  time.  Thus,  statistical  methods 
based  on  the  probabilistic  behavior  of  the  refractive  index  fluctuations  must  be  used 
to  model  the  effects  of  the  turbulence. 

Much  has  been  written  on  modeling  the  effects  of  optical  wave  propagation 
through  atmospheric  turbulence.  However,  as  Goodman  states,  “...  without  a  doubt 
the  most  influential  work  on  the  subject  is  that  of  V.I.  Tatarski”  (18:385).  This  in¬ 
fluence  is  readily  apparent  from  all  the  recent  literature  on  the  subject  that  reference 
Tatarski’s  work  (2,  3,  6,  7,  8,  9, 10, 11, 13, 15, 18,  24,  27,  37).  In  his  book,  Tatarski 
assumed  a  Kolmogorov  spectrum  (discussed  shortly)  for  $n(/c)  and  then  derived  a 
structure  function  for  the  refractive  index  fluctuations,  Dn,  as  well  as  the  amplitude 
and  phase  fluctuation’s  probability  density  functions  (PDFs)  (35). 

The  spectrum,  $n(/c),  is  a  measure  of  the  number  of  eddies  of  different  sizes 
making  up  the  turbulence  where  the  independent  variable,  ac,  is  the  wave  number. 
The  wave  number,  «,  is  inversely  related  to  eddy  size,  A,  by  /c  =  2tc/L  (recall  Figure 
1  which  illustrates  the  different  eddy  sizes)  (18:388).  Based  on  Kolmogorov’s  classic 
work  on  turbulence,  as  described  by  Coodman,  4>n(/c)  is  believed  to  consist  of  three 
regions  defined  by  the  inner  and  outer  scales  of  turbulence.  Figure  2  shows  the 
Kolmogorov  spectrum  and  its  relationship  to  the  inner,  «m,  and  outer,  k0,  scales 
and  $n(/c).  Refering  to  Figure  2,  $n(/c)  has  units  of  m3,  C2  has  units  of  m-2^3  and  n 
has  units  of  m-1.  Also  note  that  the  von  Karman  spectrum  is  often  used  to  overcome 
the  known  defects  in  the  $n(/c)  model  for  very  small  wavenumbers  (18:388-390). 

According  to  Goodman,  the  three  regions  of  turbulence  can  be  described  as 
follows.  The  first  region  consists  of  the  very  large  eddy  (or  refractive  index)  inho¬ 
mogeneities  (k  <  k0)  before  they  are  broken  up  by  the  turbulent  wind  flow  and 
convection.  The  second  region  is  the  inertial  subrange  of  the  spectrum.  This  re¬ 
gion  is  bounded  by  the  outer  and  inner  scales  of  the  turbulence  (k0  <  k  <  Km). 
It  is  in  the  inertial  subrange  that  “...  the  shape  of  the  PSD  is  determined  by  the 
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Figure  2.  Power  Spectral  Density  of  Refractive  Index  Fluctuations  (18:390) 
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physical  laws  that  govern  the  breakup  of  large  turbulent  eddies  into  smaller  ones” 
(18:389).  Finally,  the  third  region  is  the  region  outside  the  inner  scale  of  the  turbu¬ 
lence  ( k  >  Km)  in  which  the  eddies  dissipate  their  energy  due  to  the  viscous  forces 
of  the  atmosphere.  In  this  region  $n(/c)  drops  rapidly  (18:388-389). 

The  structure  function,  Dn,  is  derived  from  $n(/c)  and  describes  the  statistical 
behavior  of  the  refractive  index  fluctuations  at  two  points  in  space  by  the  equation 
(18:391) 

£>rt  =  E{[n1(n1)-n1(r^)]2}  (1) 

where 

rt-  =  a  point  in  space 

n i  =  random  fluctuations  of  the  refractive  index 
E  =  expected  value  operator 

(The  dependence  of  Dn{r^r^)  on  $„(«;)  follows  from  the  fact  that  $n(«)  is  the  three- 
dimensional  PSD  of  ni.)  Goodman,  by  assuming  homogeneous  turbulence,  solves 
Equation  (1)  to  yield  (18:392) 

D.  =  Cl  |r|2/3  (2) 

where  r  =  rx  —  r2. 

The  structure  function,  Dn,  determines  the  wave  structure  function,  D.  (The 
wave  structure  function  describes  the  statistics  of  the  optical  fluctuations  which 
result  from  propagation  through  the  randomly  inhomogeneous  turbulence.)  The 
wave  structure  function,  D ,  in  turn,  describes  the  behavior  of  the  propagating  optical 
wave  at  two  points  in  the  turbulent  medium.  The  relationship  between  Dn  and  D 
is  defined  by  (18:426) 

D  =  2.91  ib2  \r\  Dn  (3) 

where  k  is  the  optical  wave  number  ( k  =  27t/A,  A  =  wavelength).  The  wave  structure 
function  can  also  be  described  in  terms  of  the  log-amplitude  structure  function,  Dx, 
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and  the  phase  structure  function,  £>s,  by  (18:404-407) 


£(Hi,l2)  =  Dx(Xi  fotefo))  +  Ds(Si{r1),S2{r2)) 


where 

Xi{ti)  —  log-amplitude  of  the  propagating  optical  wave  at  a  point  in  space 
Si(r-)  =  phase  of  the  propagating  optical  wave  at  a  point  in  space 

■DxfxifcO.xjfe))  =  El(xita)-xita))2] 
i>s(s,(a),s2(a))  =  B[(s,(n,)  -  s2k))2] 

(Note  that,  when  considering  coherent  detection  systems,  it  is  the  phase  structure 
function,  Ds,  that  is  of  primary  interest.) 

The  wave  structure  function  is  also  related  to  the  MCF  of  an  optical  wave 
propagating  through  turbulence.  The  MCF  is  defined  as  the  cross- correlation  of  the 
optical  wave  at  two  points  in  the  observation  plane.  According  to  Lutomirski  and 
Yura,  the  MCF  describes  the  loss  of  coherence  of  a  coherently  transmitted  wave 
propagating  through  turbulence  (29:482).  In  addition,  Yura  stated  that  the  MCF 
determines  the  SNR  of  an  optical  communication  system  that  uses  coherent  detection 
(37:1399).  Lutomirski  and  Yura  concluded  by  stating  that,  under  the  conditions  of 
atmospheric  turbulence,  the  MCF  is  directly  related  to  the  wave  structure  function 
by  (29:482) 

MCF  =  exp(— Z)/2) 

The  final  quantity  of  interest  is  the  PDF  of  the  complex  amplitude  fluctuations, 
and  in  particular,  the  phase  fluctuations  of  the  propagating  optical  wave  at  a  point  in 
space.  In  his  work,  Tatarski  characterized  the  statistics  of  the  amplitude  fluctuations 
using  a  log-normal  PDF  and  the  phase  fluctuations  using  a  Gaussian  PDF  (35:Ch 
6,7). 
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With  respect  to  the  amplitude  fluctuations,  current  literature  refines  and  limits 
Tatarski’s  work  (6,  9,  7,  29).  For  example,  Churnside  and  Clifford  predict  a  log¬ 
normal  Rician  PDF  (6),  as  do  Churnside  and  Frehlich  (7).  However,  since  this 
analysis  is  concerned  with  coherent  detection,  the  amplitude  fluctuations  do  not  play 
a  major  role  in  modeling  the  optical  communication  system.  Rather,  it  is  the  phase 
fluctuations  that  dominate  the  optical  communication  system  model.  With  respect 
to  the  phase  fluctuations,  the  reviewed  literature  agrees  with  Tatarski’s  results  of  a 
Gaussian  PDF  (2,  9, 11,  13). 

2.3  Optical  Communication  Detection  Systems 

Optical  communication  detection  systems  fall  into  two  categories:  incoherent 
or  direct  detection  receivers,  and  coherent  (including  heterodyne  and  homodyne) 
detection  receivers.  Direct  detection  receivers  respond  to  the  intensity  of  the  optical 
wave,  while  coherent  receivers  respond  to  the  amplitude  and  phase. 

2.3.1  Direct  Detection  Systems  In  a  direct  detection  system,  the  signal 
current  is  proportional  to  the  received  intensity  at  the  receiver  aperture  (2,  26). 
The  direct  detection  operation  is  based  on  “photon  counting”  (25:39).  Its  primary 
advantage  is  the  simplicity  of  design  and  operation. 

2.3.2  Coherent  Detection  Systems  Coherent  receiver  systems,  on  the 
other  hand,  are  more  complex  in  design  and  operation.  In  coherent  systems,  a 
locally  generated  optical  wave  is  added  (or  mixed)  with  the  received  optical  wave 
before  photodetection  (14, 15).  This  locally  generated  optical  wave  is  equivalent  to  a 
local  oscillator  generated  waveform  in  RF  operation.  However,  there  is  a  distinction 
made  between  RF  and  optical  coherent  detection.  In  RF.  the  difference  between 
coherent  and  incoherent  detection  of  the  intermediate  frequency  (IF)  signal  depends 
on  whether  or  not  the  phase  of  the  signal  is  used.  In  optical  systems,  coherent 
detection  is  defined  as  any  detection  process  that  uses  a  local  oscillator  (26:36). 
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2.3.3  Coherent  Versus  Direct  Detection  Systems  Coherent  detection 
systems  have  several  advantages  over  direct  detection  systems.  First,  coherent 
detection  systems  offer  a  significant  improvement  in  receiver  sensitivity  (or  SNR) 
(4,  25,  26,  32).  The  SNR  improvements  stated  in  the  literature  range  from  10  dB 
(25,  26)  to  (theoretically)  20  dB  (26,  32).  Second,  coherent  detection  systems  offer 
greater  frequency  selectivity  for  the  transmitters  and  receivers  (25,  26,  32).  Third, 
they  provide  for  conveniently  tunable  receivers  (25:1).  And  finally,  coherent  detec¬ 
tion  allows  for  the  use  of  alternative  modulation  formats  (since  they  respond  to  both 
the  phase  and  amplitude  of  the  received  optical  wave)  (25:1).  However,  both  direct 
detection  and  coherent  detection  systems  are  affected  by  propagation  of  the  received 
wave  through  the  atmosphere. 

2.3.4  Atmospheric  Effects  on  Optical  Communications  Detection 
Systems  Both  direct  detection  and  coherent  detection  systems  are  adversely  af¬ 
fected  by  the  optical  wave  amplitude  and  phase  fluctuations  caused  by  atmospheric 
turbulence.  The  random  amplitude  fluctuations  result  in  reduced  gain  and  fading 
in  both  the  direct  detection  and  coherent  detection  systems  (16:13),  as  well  as  an 
overall  signal  power  loss  due  to  propagation  through  the  turbulence  (3:33). 

However,  the  coherent  detection  system  is  most  affected  by  the  random  phase 
fluctuations  of  the  optical  wave.  The  areas  of  the  optical  wave  that  are  out  of  phase 
with  the  detector  local  oscillator  will  add  destructively  resulting  in  a  reduced  IF 
amplitude  (that  is,  after  mixing)  and  thus  reduce  the  SNR  at  the  IF  stage  of  the 
detector  (16:19).  Stated  another  way,  since  the  turbulence  reduces  the  coherence  of 
the  optical  wave  (as  described  by  the  MCF),  the  effective  area  of  the  receiver  aperture 
is  reduced  to  the  size  of  the  phase  coherence  length  (11, 16).  According  to  Fried,  when 
assuming  an  infinite,  uniform  amplitude  plane  wave,  the  phase  coherence  length  is 
the  minimum  aperture  diameter  which  gives  (approximately)  the  maximum  possible 
SNR  (11:62).  Fried  (11)  and  Chen  (5)  showed  that,  by  increasing  the  aperture 
size  beyond  the  phase  coherence  length,  the  SNR  will  only  increase  by  about  3  dB. 
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Figure  3.  Normalized  Signal-to-Noise  Ratio  Versus  Normalized  Aperture  Size 

That  is,  as  Chen  states,  when  an  aperture  diameter  is  equal  to  the  phase  coherence 
length,  the  SNR  is  only  3  dB  below  what  can  be  achieved  as  the  aperture  diameter 
approaches  infinity  (5:1017).  The  relationship  between  the  phase  coherence  length, 
denoted  p0,  and  a  circular  aperture  diameter,  denoted  D,  was  shown  graphically  by 
both  Fried  (11:62)  and  Chen  (5:1017).  Chen’s  particular  representation  is  shown  in 
Figure  3. 

Lutomirski  and  Yura,  however,  found  that  Fried’s  results  for  a  maximum  useful 
receiver  diameter  are  not  valid  for  all  ranges  (29:483).  While  also  , assuming  an 
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infinite,  uniform  amplitude  plane  wave,  they  showed  that  the  SNR  increases  with 
increasing  receiver  aperture  sizes  for  path  lengths  small  when  compared  to  the  mean 
field-decay  length,  zc  (29:487).  However,  when  the  path  lengths  are  greater  than 
ze  (where  zc  is  typically  less  than  10  &m),  Lutomirski  and  Yura  showed  that  their 
results  agree  with  Pried  (29:487). 

Lutomirski  and  Buser  extended  Fried’s  results  to  the  case  of  a  finite  transmit¬ 
ting  aperture.  Their  results  can  be  summarized  as  follows: 


1.  The  largest  useful  receiver  diameter  is  the  smaller  of  ( z/k )J/2  and  the 
transverse  coherence  length,  p0,  of  the  field  at  the  receiver  due  to  the 
field  £7,4  (r)  on  the  transmitter  plane. 

2.  The  largest  useful  transmitter  diameter  is  the  smaller  of  (z/k)1/2  and 
the  transverse  coherence  length,  p0,  of  the  field  at  the  transmitter  due  to 
the  field  F(p)  in  the  receiver  plane.  (28:2159) 

The  transverse  (phase)  coherence  lengths,  p0  and  p0)  are  defined  as  the  transverse 
separations  for  which  the  correlation  of  the  propagating  optical  wave  at  two  points 
in  space  is  reduced  by  e”1,  and  where  the  wave  is  emitted  from  the  transmitter 
and  receiver,  respectively.  The  term  (z/k)1/2  defines  the  Fresnel  zone,  where  z  is  the 
propagation  path  length  (m)  and  k  is  the  optical  wave  number  (m-1)  (28:2158-2159). 

The  implications  of  Fried’s,  as  well  as  Lutomirski  and  Buser’s,  results  are  that, 
when  the  aperture  diameters  are  larger  than  both  the  Fresnel  zone  and  the  applicable 
transverse  coherence  length,  the  SNR  saturates  due  to  turbulence  (28:2159). 

It  is  interesting  to  note  here  the  varied  nomenclature  and  mathematical  ex¬ 
pressions  used  to  describe  the  phase  coherence  length.  For  example,  Fried  calls  it 
the  “efficiency  saturation  dimension”  and  denotes  it  r0  (11:57),  while  Lutomirski 
and  Buser  refer  to  it  as  the  “turbulence-induced  coherence  length”  or  “transverse 
coherence  length”  and  denote  it  p0  (28:2153,2159).  Appendix  A  presents  a  brief 
discussion  concerning  the  various  nomenclature  and  mathematical  expressions  found 
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in  the  reviewed  literature  to  describe  the  phase  coherence  length.  For  this  study, 
Yura’s  expression  is  used  (37:1403): 


p0  =  (0.545fc2C^)~‘ 


(4) 


where 

k  =  optical  wave  number  (m“x) 

C2  =  structure  constant  of  the  refractive  index  fluctuations  (turbulence  strength) 

(m~2/3) 

z  —  propagation  path  length  (m) 

(Note  that  the  terms  will  be  addressed  in  more  detail  in  Chapter  III.) 

2.3.5  Noncoherently  Combined  Large  Aperture  Optical  Heterodyne 
Receiver  Chen  studied  the  use  of  combining  multiple  receivers  as  a  means  of  over¬ 
coming  the  phase  coherence  length  limit  of  the  receiving  aperture.  His  analysis 
addressed  noncoherently  combining  the  output  signals  from  an  array  of  smaller  re¬ 
ceivers,  and  comparing  the  results  to  a  monolithic  heterodyne  receiver  of  equivalent 
aperture  size.  He  showed  that,  in  the  absence  of  turbulence,  the  monolithic  receiver 
performed  better  than  the  noncoherently  combined  small  receivers  array  when  viewed 
from  a  SNR  aspect.  However,  once  turbulence  was  introduced  into  the  analysis,  Chen 
found  that  the  noncoherently  combined  small  receivers  array  performed  significantly 
better  than  the  monolithic  receiver.  The  monolithic  receiver  was  now  severely  lim¬ 
ited  by  the  phase  coherence  length.  Chen  concluded  by  stating  that  the  optimal 
situation  occurred  when  the  receiver  diameter  of  each  of  the  small  receivers  in  the 
combined  aperture  was  on  the  order  of  the  phase  coherence  length  (5:1013). 
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2.4  Mutual  Coherence  Function  and  Signal-to-Noise  Ratio  Development 

As  stated  previously,  the  MCF  determines  the  SNR  of  an  optical  heterodyne 
communication  system.  Thus,  much  research  has  been  conducted  investigating  the 
MCF  in  the  presence  of  turbulence. 

Lutomirski  and  Yura,  by  extending  the  Huygens-Fresnel  principle,  derived  an 
expression  for  the  MCF  of  an  optical  wave  propagating  in  an  inhomogeneous  medium, 
such  as  atmospheric  turbulence  (29:1653).  According  to  Lutomirski  and  Yura 


The  principle  as  applied  to  propagation  in  a  vacuum  states  that  every 
point  of  a  wavefront  may  be  considered  as  a  center  of  a  secondary  dis¬ 
turbance  which  gives  rise  to  spherical  wavelets,  and  the  wavefront  at  a 
later  instant  may  be  regarded  as  the  envelope  of  these  wavelets... for  a 
refractive  medium,  the  extended  principle  is  that  the  secondary  wave- 
front  will  again  be  determined  by  the  envelope  of  spherical  wavelets  from 
the  primary  wavefront,  but  each  wavelet  will  now  be  determined  by  the 
propagation  of  a  spherical  wave  through  the  refractive  medium.  (29:1653) 


First,  Lutomirski  and  Yura  proved  that  the  Huygens-Fresnel  principle  could  be 
extended  to  optical  wave  propagation  in  a  weakly  inhomogeneous  medium,  such  as 
atmospheric  turbulence.  FVom  the  extended  Huygens-Fresnel  principle,  they  showed 
that  the  field  due  to  a  disturbance  across  an  aperture  can  be  calculated.  By  ap¬ 
plying  the  reciprocity  theorem  of  Helmholtz  to  the  field  and  averaging,  the  MCF 
of  the  propagating  optical  wave  results  (29:1653-1654).  The  reciprocity  theorem  of 
Helmholtz  states  that,  given  two  points  in  sp.  ce,  P0  and  Pj,  a  point  source  located  at 
P0  will  produce  the  same  effect  on  Pi  as  an  identical  point  source  at  P\  will  produce 
on  P0.  That  is,  symmetry  exists  with  respect  to  the  source  plane  and  the  observation 
plane  (17:41). 

Yura  used  Lutomirski  and  Yura’s  extended  Huygens-Fresnel  principle  and  ap¬ 
plied  it  to  a  finite  optical  beam.  Yura  showed  that  the  resulting  finite  beam  MCF  is 
valid  for  an  arbitrary  complex  disturbance  in  the  transmitting  aperture  (37:1400). 
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Lutomirski  and  Buser,  like  Yura,  also  applied  the  extended  Huygens- Fresnel 
principle  to  a  finite  optical  beam.  However,  Lutomirski  and  Buser  then  went  further 
to  apply  the  resulting  finite  beam  MCF  to  coherent  detection  (28).  First,  they  showed 
that  “a  necessary  and  sufficient  condition  for  instantaneous  reciprocity  is  that  the 
transmitter  aperture  field  and  the  receiver  weighting  be  identical”  (28:2157).  (Lu- 
tomirski  and  Buser  define  instantaneous  to  mean  for  each  realization  of  the  refractive 
index  field  (28:2157).)  They  presented  a  special  case  that  fulfills  this  condition:  a 
uniform  plane  wave  in  a  circular  aperture  transmitted  to  a  uniformly  weighted,  cir¬ 
cular  receiver  when  the  transmitting  and  receiving  apertures  are  of  equal  diameter 
(28:2157).  (The  importance  of  reciprocity  will  become  apparent  when  the  MCF  is 
derived  in  Chapter  III.) 

Lutomirski  and  Buser  then  developed  the  normalized,  average  SNR  for  coherent 
detection.  Starting  with  the  signal  power  as  developed  by  Fried  (11),  Lutomirski 
and  Buser  introduced  the  cross-correlation  of  the  complex  fields  at  two  points  in 
the  observation  plane,  px  and  P2,  due  spherical  wave  sources  at  two  points  in 
the  transmitting  aperture  plane,  ra  and  r2.  They  used  a  geometric  optics  approach 
to  derive  an  expression  for  the  cross-correlation  function,  H(p,r,z ),  where  z  is  the 
propagation  path  length,  and  p  and  r  again  represent  points  in  the  observation  plane 
and  transmitting  aperture,  respectively.  Note  that,  for  isotropic  and  homogeneous 
turbulence,  p  =  px  ~P2  and  r  =  ra  —  r2  (28:2154).  Lutomirski  and  Buser  showed  that 
the  propagation  properties  required  to  determine  the  MCF  (for  a  finite  transmitter 
source)  are  contained  in  H(p,r,z).  The  derived  signal  power  is  then  averaged,  and 
the  normalized,  average  SNR  expression  is  developed  (28:2154-2159). 

Assuming  infinite  Gaussian  apertures,  Lutomirski  and  Buser  solved  the  nor¬ 
malized  SNR  to  a  get  closed  form  solution  for  two  cases:  propagation  in  a  vacuum 
( H(p,r,z )  =  1)  and  propagation  in  turbulence  (28:2158-2159). 

The  derived  MCF  and  SNR  expressions  of  Lutomirski  and  Yura,  Yura,  Fried, 
and  Lutomirski  and  Buser  are  presented  and  discussed  further  in  Chapter  III. 
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Table  2.  Direct  Detection  Analysis  Communication  Link  Parameters  (24:57) 


Altitude(m) 

Aircraft  Speed  (m/s) 

Data  Rate  (bps) 

3048 

218.639 

20000 

3048 

218.639 

40000 

3048 

244.361 

20000 

6096 

218.639 

20000 

6069 

218.639 

IBB 

6069 

244.361 

■91 

9144 

218.639 

20000 

9144 

218.639 

40000 

9144 

244.361 

20000 

2.5  Direct  Detection  System  Bit  Error  Rate  Development 

As  discussed  in  Section  2.3.3,  coherent  detection  systems  offer  SNR  improve¬ 
ments  of  10  to  20  dB  over  direct  detection  systems.  However,  in  the  presence  of 
atmospheric  turbulence,  coherent  detection  systems  are  affected  by  both  the  ampli¬ 
tude  and  phase  fluctuations  of  the  turbulence,  while  direct  detection  systems  are 
primarily  affected  by  the  amplitude  fluctuations.  Thus,  the  SNR  improvements  co¬ 
herent  detection  offers  over  direct  detection  may  be  less  in  the  presence  of  turbulence. 

Kanavos  calculated  the  BERs  for  a  direct  detection  (noncoherent)  air-to-air 
optical  communication  link  for  various  path  lengths,  altitudes,  and  data  rates.  He 
calculated  maximum  reliable  path  lengths  for  optical  links  both  with  and  without 
turbulence  for  a  BER  of  10-5  (24:58,63).  Using  the  parameters  listed  in  Tables  2,  3 
and  4,  Kanavos  calculated  the  maximum  reliable  path  lengths  summarized  in  Table  5. 


It  should  be  noted  here  that,  although  Kanavos  does  present  his  results  in 
a  manner  similar  to  this  research  effort,  a  direct  comparison  cannot  be  made  be¬ 
tween  his  results  and  this  effort’s  results.  This  is  because  Kanavos,  in  his  research 
methodology,  requires  parameters  not  addressed  in  this  effort.  This  research,  like¬ 
wise,  requires  parameters  not  addressed  by  Kanavos. 
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Table  3.  Direct  Detection  Analysis  Transmitter  and  Receiver  Characteristics 
(24:57) 


Transmitter  Output  Power: 

50  W 

Transmitter  Optics  Efficiency: 

0.9 

Transmitter  and  Receiver  Aperture  Diameters: 

0.127  m 

Receiver  Optics  Efficiency: 

0.8 

Receiver  Gain: 

40 

Table  4.  Direct  Detection  C%  Summary  (24:77) 


Altitude  (m) 

Cl  (m-!/3) 

3048 

6.4  xlO"17 

6096 

1.617  x  10"17 

9144 

1.947  x  lO"18 

Table  5.  Direct  Detection  Analysis  Results  Summary  (24:73) 


Maximum  Reliable  Path  Lengths 

Altitude  (m) 

Baseline  Link  (km) 

Turbulence  Link  (km) 

3048 

95 

42 

6096 

208 

77 

9144 

239 

134 
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2.6  Summary 

The  literature  indicates  that  optical  communication  is  possible  through  the 
atmosphere.  However,  the  atmosphere  does  degrade  the  performance  of  the  optical 
communication  system. 

Tatarski’s  work  in  wave  propagation  through  turbulence  laid  the  foundation 
for  determining  the  effects  of  atmospheric  turbulence  on  an  optical  wave  propagating 
through  the  atmosphere.  Subsequent  literature  refined  Tatarski’s  results,  particu¬ 
larly  in  the  area  of  the  amplitude  fluctuation  PDF  of  the  optical  wave. 

The  effects  of  the  turbulence  also  affect  the  detection  of  the  optical  wave. 
Without  turbulence,  the  literature  agrees  that  coherent  detection  offers  better  re¬ 
ceiver  sensitivity  (or  SNR)  than  direct  detection.  However,  although  both  coherent 
and  direct  detection  systems  are  affected  by  the  random  amplitude  fluctuations  of 
the  optical  wave,  the  coherent  detection  system  is  further  affected  by  the  phase 
fluctuations. 

Bit  error  rate  calculations  for  direct  detection  systems  indicate  path  lengths 
of  approximately  40  to  134  km  are  possible  (for  a  BER  of  10-5)  in  the  presence  of 
turbulence. 
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III.  Theory  and  Methodology 


This  chapter  addresses  three  topics:  general  SNR  equation  development,  op¬ 
timum  aperture  dimensions,  and  uniformly  weighted  aperture  SNR  and  BER  equa¬ 
tions  development.  Specifically,  this  chapter  presents  the  development  of  the  SNR 
equation  for  coherent  optical  communications  in  the  presence  of  turbulence.  Second, 
a  closed-form  solution  for  the  SNR  is  obtained  by  making  assumptions  about  the 
aperture  weighting  functions  (Gaussian)  and  cross-correlation  function  (Gaussian). 
From  this  closed-form  solution,  an  analysis  of  SNR  versus  aperture  dimensions  is 
made.  Third,  making  new  assumptions  about  the  aperture  weighting  functions  (uni¬ 
form),  an  integral  expression  for  the  SNR  is  developed.  This  integral  expression  is 
solved  using  numerical  integration  methods.  This  chapter  concludes  with  a  short 
summary. 

3.1  General  Signal-to-Noise  Ratio  Development 

In  this  section,  the  general  SNR  for  coherent  optical  communications  in  the 
presence  of  turbulence  is  developed.  First,  the  Huygens- Fresnel  principle  is  presented 
for  vacuum  propagation.  Using  the  work  of  Lutomirski  and  Yura,  the  Huygens- 
Fresnel  principle  is  then  extended  to  the  application  of  propagation  in  a  turbulent 
medium.  Second,  the  MCF  for  optical  propagation  in  a  turbulent  medium  is  devel¬ 
oped.  This  development  includes  the  development  of  the  cross-correlation  function, 
H(p,r,  z).  Finally,  the  general  SNR  equation  is  developed. 

3.1.1  Huygens-Fresnel  Principle  Before  discussing  the  Huygens- Fresnel 
principle  for  vacuum  propagation,  a  brief  review  of  optical  propagation  for  monochro¬ 
matic  light  is  presented.  First,  let  the  function  u(P,t)  represent  “...the  scalar  ampli¬ 
tude  of  one  polarization  component  of  the  electric  or  magnetic  field  associated  with 
a  monochromatic  optical  disturbance”  where  P  is  some  point  in  space  and  t  is  some 
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instant  in  time  (18:117).  An  analytic  signal  representation  associated  with  u(P,t) 
can  be  written  as  (18:117) 

u(P,  t)  =  U[P)  exp  (—ikct)  (5) 


where 

£7(P)  =  phasor  amplitude 
k  =  optical  wave  number  (m) 
c  =  speed  of  light  =  3  x  108  m/s 
t  =  time  (s) 

Referring  to  Figure  4,  the  effect  of  the  wave,  u(P,i),  on  an  infinite  surface  E, 
as  well  as  the  effect  of  the  wave  at  some  point  beyond  E,  can  be  examined.  Let 
u(P,t)  be  incident  on  E  from  the  left.  To  determine  the  phasor  amplitude,  f/(P),  at 
some  point  to  the  right  of  E,  say  P0,  in  terms  of  the  field  on  E,  the  Huygens-Fresnel 
principle  can  be  applied.  Note  that,  in  Figure  4,  6  is  the  angle  between  the  normal  to 
E  and  the  line  joining  P\  and  P0,  and  z  is  the  distance  between  Pi  and  P0  (18:117). 

The  Huygens-Fresnel  principle  is  based  on  the  theory  that,  given  a  spherical 
wavefront  of  an  optical  disturbance,  if  each  point  on  that  wavefront  were  taken  to 
be  a  “secondary”  spherical  disturbance,  then  the  total  wavefront  at  any  later  instant 
in  time  could  be  constructed  from  the  envelope  of  the  “secondary”  wavelets  (17:31). 
Mathematically  (and  refering  to  Figure  4),  the  Huygens-Fresnel  principle  can  be 
expressed,  for  z  A,  as  (18:117) 

(6) 


where 

A  =  wavelength  (m) 
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5  =  infinite  surface  on  which  the  wave  is  incident 
z  —  distance  between  P\  and  P0  (m) 

6  =  angle  between  the  normal  to  S  and  the  line  joining  P0  and  P\ 

X  (0)  =  obliquity  factor 

dS  =  elemental  surface  area  on  S 

Lutomirski  and  Yura  extended  the  Huygens- Fresnel  principle  for  vacuum  prop¬ 
agation  to  optical  propagation  in  a  weakly  inhomogeneous  medium  (such  as  atmo¬ 
spheric  turbulence)  (29:1653).  In  addition,  they  considered  the  disturbance  as  spec¬ 
ified  over  some  finite  aperture,  as  opposed  to  infinite  plane  waves,  infinite  spherical 
waves  or  infinite  Gaussian  disturbances  (29:1652). 

Under  the  consideration  of  a  finite  aperture,  a  new  parameter,  r,  is  introduced 
to  denote  a  point  in  the  transmitter  plane,  while  the  parameter  p  will  now  denote 
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a  point  in  the  observation,  or  receiver,  plane.  The  transmitter  and  receiver  planes 
are  assumed  to  be  parallel  and  lie  in  an  x-y  orientation  (where  z  is  the  axis  of 
propagation). 

Assuming  the  scattering  angles  are  sufficiently  small  (that  is,  0  «  0  so  that 
x(9)  ~  1),  the  aperture  lies  in  a  plane  normal  to  the  z  axis  (axis  of  propagation), 
and  the  quantity  k  |p  —  r|  1,  Lutomirski  and  Yura  show  (29:1653): 
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where 


ip  (p,  rj  describes  the  turbulence-induced  phase  perturbations  on  a  spherical  wave 
propagating  from  a  point  r  to  a  point  £  (rad)(37:1401) 

es  =  unit  vector  from  r  to  p 

en  =  unit  vector  in  direction  normal  to  the  wavefront  in  the  transmitter  plane 

e2  =  unit  vector  in  direction  along  the  z  axis 

d2r  =  elemental  area  at  the  point  r  in  the  aperture 


When  the  distance  of  propagation  is  large  with  respect  to  the  aperture  size 
(paraxial  approximation),  then  U_  (p)  reduces  to  (29:1654) 


(8) 


where 


z  =  distance  between  the  points  r  (in  transmitter  plane)  and  p  (in  receiver  plane) 


z0  =  perpendicular  distance  between  the  transmitter  and  receiver  planes 
Oq  =  angle  between  the  normal  to  the  transmitter  plane  and  r  —  p 


The  geometry  is  illustrated  in  Figure  5. 

Equation  (8)  can  be  reduced  and  rewritten  using  the  following  definitions  and 
assumptions  (37:1400): 


1.  Define  Ua(l)  to  be  the  arbitrary  complex  disturbance  at  r,  where  the 
subscript  A  implies  the  disturbance  defined  in  the  transmitter  plane. 

2.  Assume  60  ~  0. 

3.  Define  the  field  at  p  due  to  a  spherical  wave  propagating  in  the  refractive 
medium  from  r  as 


exp 
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(9) 


Note  that  G  (p,r)  is  simply  one  representation  of  a  Green’s  function. 
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Now  Equation  (8)  becomes  (37:1400) 

V[p)=~jG(p,r)UA{z)<Pr  (10) 

where  the  integration  is  carried  out  over  the  transmitter  aperture  (37:1400). 

8.1,2  Mutual  Coherence  Function  Development  In  this  section,  the 
field  amplitude  equation  developed  in  Section  3.1.1  is  used  to  determine  the  MCF 
of  the  optical  waves  in  the  presence  of  atmospheric  turbulence.  The  MCF  is  then 
manipulated  to  “pull  out”  the  atmospheric  effects  in  terms  of  a  cross-correlation 
function  H(p>  r,  z).  Recall  that  H(p,r ,  z)  was  defined  to  be  the  two  point  correlation 
function  (28:2154).  The  relationship  between  H(p,r,z )  and  the  phase  structure 
function,  Ds ,  is  then  presented. 

The  MCF,  which  represents  a  cross-correlation  of  the  propagating  optical  wave 
at  points  p ^  and  p2  is  defined  to  be  (18:174) 

MCF^{v(pt)u:(p2))  (11) 

where  <  •  >  represents  the  ensemble  average.  Assuming  U (p.)  is  the  field  due  to  a 
transmitting  aperture  with  field  distribution  Ua(j^),  Equation  10  is  substituted  into 
Equation  11  to  yield 

MCF  =  (^)  /  /  (G  (Hi’-1)  Ua  ta) G*  U\  (r2))  dVjdVj  (12) 

Note  that  <  Ua  (zj)  U\  (£2)  >=  Ua  (£1)  Ua  (£2)*  thus  the  quantity  can  come  out  of 
the  ensemble  average,  which  yields  the  MCF  derived  by  Yura  (37:1400): 

MCF  =  (^r)  / /  (G  UA  (Li)  Ua  {ra)#rxd2r2  (13) 
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Lutomirski  and  Yura  proved  that  the  reciprocity  theorem  holds  for  G  (g,  r) 
(29:1658),  that  is  the  field  at  p  due  to  a  point  source  at  r  is  identical  to  the  field  at 
r  due  to  a  point  source  at  p,  thus  (37:1401) 

a  fo.r,)  G-  (&,r2)  =  G  (r„&)  G‘  (*,£,)  (14) 


In  developing  the  expression  for  the  MCF  as  a  function  of  a  cross-correlation  func¬ 
tion,  H(p,r,z )  (as  defined  by  Yura  (37:1401)),  the  importance  of  the  role  that  the 
reciprocity  theorem  plays  becomes  evident.  Recalling  Equation  (9),  then  (37:1401) 
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where 

H  (pvLdP2,!-2>z)  -  (exp  [rp  (ei.pJ  +  (l2,P2)])  (16) 

Thus  adherence  by  the  MCF  to  the  reciprocity  theorem  allows  the  phase  perturba¬ 
tions,  V’foiE,)?  to  be  expressed  in  a  form  consistent  with  Yura’s  (37:1401)  definition 
of  H  z)-  (That  is,  if  reciprocity  did  not  hold,  ^(r,-,^)  would  not  neces¬ 

sarily  equal  ^(p^r,).  The  cross-correlation  function,  H  (Pj > Hi > g2 > 2I2 >  ^) >  35  defined 
by  Yura,  could  not  then  be  used  in  the  MCF  defined  here.) 

Recalling  the  assumption  that  z ,  the  distance  of  propagation,  is  large  compared 
to  the  aperture  size  (known  as  the  paraxial  approximation),  then  the  |rj  —  pj  and 
1-2  -  £2 1  terms  in  the  denominator  of  Equation  (15)  can  be  replaced  by  z  to  yield 


(G(Li,Pi)G-(r2,p2)) 


^exp  [ik  (|ei  —  £j| 
xif  (£1,r1,p2,r2,z) 


(17) 
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Yura  shows  that,  again  using  the  paraxial  approximation  (37:1401), 
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(18) 


which  now  yields 


(0(r,,Ej0-  (&,&))  = 


x^(£1>Ei>£2>E2^) 
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(19) 


Now  an  expression  for  the  MCF  can  be  derived  by  substituting  Equation  (19) 
into  (13)  to  yield 


MCF  = 


J  Jexp{rzh~r-' 


.  2nz . 


V2-ta 


(g,, Ua  (u)U\  (r2)rp n<P& 


(20) 


where  the  terms  are  previously  defined.  For  isotropic  and  homogeneous  turbulence, 
H  (px, £1,^5, £2,  is  only  a  function  of  the  difference  coordinates  so  that  (28:2154): 


mcf  =  (2^)  //exp{s[i2<-&r-b-&r]} 

xH[h~  &>&  “  E».*)  Kl  (&)  (is) d1  Ltd  Li 


(21) 


Now  that  a  general  expression  for  the  MCF  has  been  derived,  the  function 
H  fpj  —  p2,li  —  r2,«)  can  be  investigated.  Recall  Equation  (16)  which  states 

H  (gj  -  p2, tx-r2,z)  =  (exp  [V>  (rx,  pj  +  (r2, p2)] )  (22) 

Yura  solved  Equation  (22)  through  second  order  in  nx  where  ^>(r,-,p.)  is  related  to  nx 
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via  the  path  integral  of  ni  along  the  optical  path  from  p.  to  r,-.  Recall  that  V’(n»>2.) 
describes  the  turbulence-induced  phase  perturbations  of  the  wave,  and  that  nj  is  a 
measure  of  the  refractive  index  fluctuations  about  the  free-space  value  (where  the 
free-space  value  is  approximately  unity).  Without  presenting  the  details  of  Yura’s 
work  here,  Yura  determined  that  the  general  form  of  H  (px  —  p2,t\  —  £2,  z'j  can  be 
expressed  as  (37:1406) 

H  (pj  - -rj ,z)  =  exp {-  [(|V>i)|2  -  (V’l  (li^)  V>i*  (e2,£2))]}  (23) 

where 

(iV’il2)  =  2tt  k2z  J  $n  (/c)  dn  (24) 

Yura  shows  that  Equation  (24)  is  correct  through  second  order  in  ni,  but  also 
points  out  that,  given  the  experimental  evidence  and  various  theoretical  arguments 
showing  ^  to  be  a  Gaussian  random  variable,  Equation  (24)  may  be  correct  to  all 
orders  in  n\  (37:1401). 

Note  that,  for  simplicity,  the  implicit  dependence  of  H  —  P2>Ii  —  r2,  z)  on 
z  is  dropped  for  the  remainder  of  this  study. 

3.1.3  Signal-to-Noise  Ratio  Dependence  In  this  section,  the  relation¬ 
ship  between  the  MCF  and  the  signal  power,  S’,  is  presented.  The  signal  power  is 
averaged  and  then  divided  by  the  noise  power,  JV,  to  obtain  a  SNR.  This  SNR  is  then 
normalized  by  the  effective  cw  transmitter  output  power,  W The  section  concludes 
with  an  expression  for  the  general  normalized  SNR. 

Lutomirski  and  Buser  show  that  the  MCF,  as  defined  by  MCF  =  (U{p1)lT(p2)\ 
is  related  to  the  signal  power,  5,  by  (28:2157) 

S  =  |  {r)ARgf  R  |  J  U(p)F  (p)  d2p |  (25) 
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where 


tj  =  quantum  efficiency  of  the  receiver  photodetector 

Ar  =  amplitude  of  the  receiver  local  oscillator  signal 

g  =  gain  of  the  receiver  current  amplifier 

R  =  resistance  of  the  receiver  load 

F(p)  =  receiver  aperture  weighting  function 

Substituting  Equation  (21)  into  (25)  and  taking  the  ensemble  average  yields 
(28:2157) 

<s>  =  ff 

xH  ~  Li)  <Pli #Lid\d2p2  (26) 

Lutomirski  and  Buser  (assuming  circular  apertures)  stated  that,  for  large  cur¬ 
rent  amplifier  gains,  g,  the  noise  affecting  the  system  will  be  dominated  by  shot  noise 
(since  Ar  is  much  larger  than  the  received  signal  amplitude),  thus  (28:2158) 

N  =  eg*i1(wH)DlAiRR  (27) 


where 

N  =  noise  power  at  the  receiver  per  unit  bandwidth  (W  •  m2/Hz) 
e  =  charge  on  an  electron  =  1.602  x  10-19  ( coul ) 

Dr  =  receiver  diameter  (m) 

Lutomirski  and  Buser  also  pointed  out  that,  for  the  problem  of  cw  lasers  (as 
in  this  analysis),  multiplying  the  noise,  N,  by  the  effective  cw  transmitter  output 
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power,  Wj 4,  will  result  in  a  final  SNR  normalized  to  1  watt  of  cw  transmitter  output 
power.  They  defined  WA  to  be  (28:2158) 

WA  =  J\UA(r)\2<Pt  (28) 

They  then  went  one  step  further  to  define  a  constant  B  such  that 

WA  =  BD\  (29) 

where 

B  —  a  constant  for  a  given  aperture  field  UA(r) 

Dt  =  transmitter  diameter  (m) 


Thus,  combining  Equations  (26),  (27)  and  (28)  yields  the  final  result  of  this 
section,  the  general  SNR  (28:2158): 


(S) 

nwa 


h  (!)  W 

x«tpg[i&-&r-b-4]} 

xH{h~  £2>-i  ”  &)  (30) 


where 

V'A  (r.)  V'x  fa)  =  j &UA  fa)  U-A  fa)  (31) 

Lutomirski  and  Buser  (28)  solved  the  general  SNR  expression,  Equation  (30), 
while  assuming  infinite  Gaussian  weighting  functions  for  circular  transmitter  and 
receiver  apertures.  Although  convenient  in  that  they  provide  for  a  closed  form  so¬ 
lution  of  the  general  SNR  equation,  the  infinite  Gaussian  weighting  functions  are 
not  realizable  due  to  their  infinite  extent.  Fried  accomplished  his  analysis  of  optical 
heterodyne  detection  in  the  presence  of  atmospheric  turbulence  while  assuming  uni- 
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formly  weighted,  circular  receiving  apertures  (11:58).  However,  Pried  assumed  an 
infinite  plane  wave  incident  on  the  receiver.  For  this  research  effort  a  finite  wave  is 
assumed  since  a  finite  transmitter  aperture  is  used.  The  finite  aperture  is  a  realizable 
scenario,  unlike  Fried’s  infinite  plane  wave  assumption. 

Unfortunately,  as  discussed  in  Appendix  B,  uniformly  weighted  circular  aper¬ 
tures,  when  used  in  the  general  SNR  equation,  Equation  (30),  result  in  a  SNR 
expression  consisting  of  a  single  four-fold  integral.  The  problem  with  the  four-fold 
integral  is  that  it  cannot  easily  be  reduced  analytically.  To  use  numerical  methods 
to  solve  the  SNR  equation,  the  SNR  must  be  reduced  to  an  expression  consisting  of 
fewer  than  four  integrals  (due  to  the  limited  processing  capability  of  the  computer’s 
available).  Thus,  in  order  to  solve  the  SNR  for  the  purposes  of  this  feasibility  study, 
square  apertures  have  been  assumed.  As  will  be  shown,  square  apertures  result  in  a 
SNR  expression  consisting  of  2  two-fold  integrals. 

The  dimensions  of  the  apertures,  to  this  point,  have  been  with  respect  to  the 
diameter  of  a  circular  aperture.  To  make  meaningful  comparisons  from  a  feasibility 
standpoint,  the  dimension  of  the  square  aperture  will  be  calculated  by  considering 
a  circular  aperture  diameter,  and  then  determining  the  associated  square  aperture 
dimension  by  equating  the  areas  of  the  uniformly  weighted  apertures. 

The  following  sections  solve  Equation  (30)  for  two  cases.  First,  Lutomirski  and 
Buser’s  approach  (infinite  Gaussian  weighting  functions  for  the  circular  transmitter 
and  receiver  apertures)  is  developed.  This  development  is  accomplished  for  two 
purposes:  (1)  it  provides  a  baseline  for  the  SNR  work  that  has  been  accomplished 
in  the  past,  and  (2)  it  provides  a  simple,  closed  form  solution  to  the  general  SNR 
expression.  This  closed  form  solution  can  then  be  evaluated  to  determine  trends 
upon  which  the  following  approach  can  be  compared. 

The  second  approach  involves  extending  the  results  of  the  general  SNR  expres¬ 
sion  to  uniformly  weighted,  finite,  square  apertures.  This  derivation,  as  discussed 
above,  provides  for  a  more  realistic  optical  communication  system  than  either  Fried’s 
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(11)  or  Lutomirski  and  Buser’s  (28)  scenarios.  In  addition,  this  square  aperture  SNR 
will  provide  the  basis  for  the  feasibility  evaluation. 

It  should  again  be  recalled  that  this  study  is  limited  in  that  it  only  addresses  the 
losses  due  to  propagation  path  length  and  turbulence-induced  phase  perturbations. 


3.2  Effective  Aperture  Sizes  Development 

In  this  section,  a  closed  form  solution  for  the  general  SNR  expression,  Equa¬ 
tion  (30),  is  derived.  From  this  SNR  expression,  the  relationship  between  the  SNR 
and  the  phase  coherence  length,  p0,  is  evaluated.  The  results  of  this  evaluation  are 
used  in  a  trends  analysis  of  the  SNR  versus  aperture  dimensions.  This  data  is  also 
used  in  a  preliminary  investigation  of  the  effect  of  wavelength  on  the  SNR.  Recall 
that  p0  describes  the  separation  distance  at  which  the  correlation  between  two  points 
of  a  propagating  optical  wave  is  reduced  by  e-1. 

As  developed  in  Section  3.1,  the  SNR  normalized  to  1  watt  of  cw  output  power 
can  be  written 
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where 

7 i  =  quantum  efficiency  of  the  photodetector  (receiver) 
e  =  charge  on  an  electron  ( coulombs ) 
k  =  wavenumber  =  y-  (m-1) 

WA=f\UA(r)\2d?r  (W  •  m2) 

Dr  =  receiver  aperture  diameter  (m) 
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z  =  propagation  path  length  (m) 

r;  •-=  vector  in  the  source  plane  (transmitter)  (m) 

p.  =  vector  in  the  observation  plane  (receiver)  (m) 

UA(r)  =  transmitter  aperture  field  distribution 
F(p)  --  receiver  aperture  distribution 

H(p i  —  P2,£i  —  r2)  =  cross-correlation  function  (as  defined  by  Equation  (23)) 

Assuming  nontruncated  Gaussian  weighting  functions  for  the  circular  trans¬ 
mitter  and  receiver  apertures,  UA(r)  and  F(p)  can  be  written 

u*l B  = 

Hit  =  exP  (2^)  (33) 

where 

U0  =  field  at  the  center  of  the  transmitting  aperture 
a  =  effective  transmitter  aperture  diameter  (m) 
b  =  effective  receiver  aperture  diameter  (m) 

This  results  in 

w*=ijy  Miy)]’*  <34) 

Using  the  fact  that  r  =  rxx  +  ryy,  and  the  definite  integral  relation  (20:307) 

/  exp  (— p2x2)  dx  =  (35) 

J- 00  v  '  p 

then 

WA  =  itU20a2  (36) 
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With  respect  to  H  {px  —  P2 ’  — 1  —  £2)’  Lutomirski  and  Buser  state  that,  when 
the  spherical  wave  coherence  function  is  Gaussian,  the  corresponding  cross-correlation 
function  is  given  by  (28:2158) 


=exp  —^5  |(p,  —  g2)  —  (ni  —  E2) 


(37) 


where  p0  =  phase  coherence  length  (m).  The  use  of  a  Gaussian  distribution  for 
H  ^  —  £2,£i  —  r2)  is  discussed  in  more  detail  in  Section  3.3.2. 

Substituting  Equations  (33),  (36)  and  (37)  into  (32)  and  integrating  the  quadru¬ 
ple  integral  over  (—00, 00)  yields  (28:2159) 
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where  C  =  8  (?) 

Equation  (38)  can  be  reduced  further  by  assuming  equal  effective  transmitter 
and  receiver  aperture  diameters  (a  =  b  =  d ).  Thus  Equation  (38)  is  reduced  to 
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where  d  =  effective  transmitter/receiver  aperture  diameter  (m). 
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3.3  Signal-to-Noise  Ratio  and  Bit  Error  Rate  Development 

In  this  section,  a  SNR  expression  for  uniformly  weighted,  finite  apertures 
is  derived.  First,  the  expressions  for  the  noise  power,  JV,  and  the  normalization 
factor,  Wj 4,  are  derived.  (Note  that  the  noise  power,  TV,  will  differ  slightly  from 
that  presented  in  Section  3.1.3  (Equation  (27))  due  to  the  square  apertures.)  Sec- 
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ond,  the  SNR  equation  for  square  apertures  is  presented.  Following  this  presenta¬ 
tion  of  the  SNR,  an  assumption  is  made  regarding  the  cross-correlation  function, 
H(px  —  P2,Li  — 12))  and  variable  substitutions  are  made  in  order  to  reduce  the  SNR 
to  an  expression  which  allows  numerical  integration  methods  to  be  used  to  solve  the 
given  expression.  After  the  SNR  equation  is  derived,  the  BER  equation  is  developed. 


3.3.1  Noise  Power  and  Normalization  Factor  Under  the  assumption  of 
uniformly  weighted,  square  apertures,  the  receiver  and  transmitter  aperture  distri¬ 
butions  can  be  written 


Ua(l) 

=  U,W(r) 

(40) 

F(t) 

=  Wfe) 

(41) 

where 


U0  =  amplitude  of  the  transmitting  aperture  field 
VF(ut  )  is  defined  as 

VF(&)  = 


* 

1 

0 

<? 


where  v  represents  either  r  or  p  and  i  represents  the  x  and  y  coordinates. 


Under  the  assumption  of  square  apertures,  minor  adjustments  must  be  made  to 
the  noise  equation,  Equation  (27)  in  Section  3.1.3,  which  assumed  circular  apertures. 
Thus,  Equation  (27)  becomes 


N  =  eg^DR=DR..A2RR 


(42) 


where 

Drx  =  receiver  aperture  length  in  the  x-direction  (m) 
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Dfty  =  receiver  aperture  length  in  the  ^/-direction  (m) 


Equation  (28)  in  Section  3.1.3  can  now  be  evaluated  as  follows: 
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where 

Dtx  =  transmitter  aperture  length  in  the  x-direction  (m) 
Dxy  =  transmitter  aperture  length  in  the  y-direction  (m) 


Since  square  apertures  are  being  assumed,  Equations  (42)  and  (43)  can  further 
be  refined  as 


N  =  eg2r]DlA2RR 

(44) 

$ 

II 

b 

(45) 

where 


Dr  —  Drx  =  Drv 

Dt  =  DTx  =  Dfy 


3.3.2  Square  Aperture  Signal-to-Noise  Ratio  Development  The  SNR 
equation  for  square  apertures  can  be  derived  by  combining  Equations  (30),  (44)  and 
(45)  to  yield 
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To  solve  this  expression,  a  closed  form  solution  for  H{px  —  p_2,L\  —  r2)  must 
first  be  determined.  Several  authors  present  solutions  for  —  £2,£i  —  21a))  hut 
these  solutions  do  not  all  agree.  For  example,  Yura  solved  the  general  expression  for 
if(p,r),  Equation  (23),  by  assuming  a  Kolmogorov  spectrum  for  $„(«).  His  resulting 
expression  is  (37:1402) 


H(Ei  “£2>£i  -£2)  =  exp< 


2.91,o^o 

—kch 


/3\  (!&-&!’- fa -&l§) 

U/  |(Ei_E2)_(!;i_r2)| 


(47) 


Barakat  and  Beletic,  on  the  other  hand,  dispute  the  use  of  the  Kolmogorov 
spectrum  for  $n(/c)  (1)-  According  to  Barakat  and  Beletic  (as  well  as  sources  they 
cite  in  their  paper),  the  Kolmogorov  approach  “...does  not  satisfy  the  necessary  dif¬ 
ferentiability  conditions  at  s  =  0  and  thus  cannot  represent  a  spatially  stationary, 
isotropic  surface”  (1:669),  where  s  is  defined  to  be  [(pi  —  p2)2  +  (ri  —  r2)2]1^2  (1:658). 
They  chose  to  approximate  the  phase  structure  function  with  a  Gaussian  distribu¬ 
tion.  To  quote  a  statement  they  attribute  to  another  paper  by  Barakat,  (noting  that 
their  variable  W3(^rj)  is  equivalent  to  Ds{p,L )) 


In  view  of  the  fact  that  no  one  has  really  been  able  to  specify  the  sta¬ 
tistical  properties  of  W3((,ri)  by  theory  or  by  experiment,  it  helps  to 
make  some  plausible  assumptions  of  the  statistical  behavior  of  W4(£,  77) 
consistent  with  the  physics  of  the  situation.  W4(£,t?)  is  taken  to  possess 
an  underlying  Gaussian  probability  structure.  This  is  purely  a  working 
hypothesis  backed  in  part  by  a  central  limit  theorem  argument.  (1:654) 

As  a  result  of  Barakat  and  Beletic’s  assertions,  as  well  as  the  fact  that  Yura’a 
expression  for  H(px  —  p^ti  —  £2)  1S  not  separable  due  to  the  j(gj  —  p2)  —  (rj  —  r2)| 
term  in  the  denominator  of  the  exponential  function  (Equation  (47)),  a  Gaussian 
cross-correlation  function  is  assumed  for  this  study.  It  should  be  noted  here  that 
the  use  of  a  Gaussian  representation  for  the  cross-correlation  function  is  common 
practice.  For  example,  Lutomirski  and  Buser,  as  shown  in  Section  3.2,  assumed  a 
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Table  6.  Five-Thirds  Versus  Gaussian  Cross-Correlation  Function  Parameters 


A(nm) 

| 

1 

z(km) 

1.5 

■9 

830 

•Bill  S3 

wM 

904 

100 

Gaussian  cross-correlation  function  when  deriving  the  closed  form  solution  to  Equa¬ 
tion  (30).  For  this  study,  assuming  a  Gaussian  cross-correlation  function  allows  the 
SNR  of  Equation  (46)  to  be  reduced  to  an  expression  consisting  of  2  two-fold  integrals 
(as  stated  earlier). 

The  Gaussian  cross-correlation  function  can  be  written  (28:2158) 


H  {h  ~  ~  L*)  =exP 


-fife -£,)-(* -4 


(48) 


where  p0  is  the  phase  coherence  length  and  is  defined  by  Equation  (4).  (See  Appendix 
A  for  a  discussion  of  the  phase  coherence  length.) 

Note  that  the  assumption  of  square  apertures  is  adequate  from  a  feasibility 
standpoint.  That  is,  in  the  analysis  of  the  feasibility  of  coherent  optical  communica¬ 
tions  in  an  air-to-air  role,  the  shape  of  the  aperture  will  not  affect  the  analysis.  Note 
also  that  the  Gaussian  cross-correlation  function  assumption  closely  approximates 
Yura’s  five-thirds  cross-correlation  function  in  most  cases.  Figures  6  through  8  illus¬ 
trate  examples  of  the  two  cross-correlation  functions  for  wavelengths  of  1.5,  830  and 
904  nm,  respectively,  with  the  specific  parameters  used  to  plot  the  functions  listed 
in  Table  6.  In  all  cases,  the  difference  r_\  —  £2  assumed  to  be  zero. 

With  the  Gaussian  cross-correlation  function,  the  SNR  equation  can  now  be 
written 
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Figure  8.  Comparisons  of  the  Five-Thirds  and  Gaussian  Cross-Correlation  Func¬ 
tions:  A  =904  nm 


<  S  >  _  k2  /r)\  1 

J  J  J  J  w  (h)W- (rjw  (£l)  w  y 

xexp  -i|(g, -£j)  -  (n, -rj)|2  &vJ?vJ?iL\£u  (49) 
.  Po 

where 

p0  =  (0.545Ai2C^)’’?  (m) 

C2  =  refractive  index  structure  constant  (m-2/3) 

Note  that  the  structure  constant,  C2,  is  generally  a  function  of  altitude  and  is  often 
expressed  analytically  using  the  Hufnagel  model  (23:6-14): 

Cl  =  8.2  x  10-"W°exp  (-^j)  +2.7  x  lCrKexp  (~^J  (50) 


0.02  0.04  0.06  6.08  0.1  0.12  0.14 

d  (cm) 
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where 


h  =  altitude  above  sea  level  (m) 

V  =  relative  rms  wind  speed  (aircraft  velocity)  (m/s) 

In  order  to  reduce  Equation  (49),  sum  and  difference  variable  substitutions  axe  made 
as  follows: 


Al;  = 


Ex,-  = 


so  that 


P{  -  Sxi  +  ^Ax,- 

u  =  au-jAasi  (52) 


The  Jacobian  for  this  transformation  is  unity. 

Substituting  Equations  (51)  and  (52)  into  (49),  and  rearranging  the  integrals 

yields 


<  S  > 
NWa 


k 2 


87 r222 


X  / /exp  (|Ai,|2  -  |Ai,|!)j  exp  j-i  [|Ai,  -  Ax/]  j 

{J]W  (Ss,  -  jAi,)  W  -  iAx2) 


x  W  (Sxt  +  jAi,)  W  (Eij  +  jAs2) 
xd2  Ax1d2Ax2 


(53) 
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Rewriting  the  vectors  in  terms  of  their  x  and  y  components 


X j  —  Axzi«r  “f"  AtTyiJ/ 

S Xi  =  S  xx{X  +  Ezyij/  (54) 

The  inner  integrals  of  Equation  (53)  can  be  solved  using  a  geometric  approach. 
Assuming  equal  sized,  square  transmitter  and  receiver  apertures,  a  new  variable,  D, 
can  be  defined  such  that  D  =  Dj  —  Dr.  Note  that  this  assumption  follows  from 
use  of  the  reciprocity  theorem.  Recall  that,  as  stated  explicitly  by  Lutomirski  and 
Buser,  “...a  necessary  and  sufficient  condition  for  instantaneous  reciprocity  is  that 
.the  transmitter  aperture  field  and  the  receiver  weighting  be  identical  (28:2157). 
That  is,  the  transmitter  and  receiver  apertures  must  be  equal  in  size  and  have  the 
same  form  of  aperture  weighting  (in  this  case,  uniform  weighting). 

Figure  9  shows  the  geometry  used  to  solve 


xW  x  +  (pXyi  +  ^Axj,^  yj 

xdT/XxidSxyi  (55) 

As  shown  in  Figure  9,  the  area  of  overlap  is 

Area  =  (D  -  |Axzt|)  ( D  -  |Axy,|)  (56) 


where 


—D  <  A xxi  <  D 

-D  <Axj „•<  D  (57) 
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Figure  9.  Sum  Integral  Geometry 


Thus,  the  SNR  now  becomes 
<S> 

nwa 


The  above  SNR  can  now  be  solved  using  numerical  methods.  Appendix  C 
discusses  the  numerical  methods  used  to  solve  the  SNR  equation  and  the  resulting 
computer  code. 

3.3.3  Bit  Error  Rate  Development  In  this  section,  the  BER  equation  will 
be  developed.  First,  the  BER  will  be  defined  for  optical  communication  systems. 


-  is?  (!)  £  C  L  <D  - 

xexp  £  (Al^i  - 

x  exp  r  f  dAx^i  dAr  ^2 

Po  J 

X  I-D  I-D  ~  exp  Yz  (A*S*  “  Ax^2) 

x  exp  \~~fi  (^xv2  -  Axj,i)2 
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Then  the  relationship  between  the  BER  and  the  SNR  will  be  discussed.  The  section 
will  conclude  with  the  BER  expression  used  for  this  analysis. 

The  BER,  for  binary  systems,  is  equivalent  to  the  probability  of  error,  P(e). 
Assuming  the  number  of  ones  and  zeros  transmitted  are  the  same,  then  the  proba¬ 
bility  of  error  becomes  (33:430) 

P(e)  =  i|P(0/l)  +  f(l/D)]  (59) 


where 


P(0/1)  =  probability  of  a  “0”  being  received  when  a  “1”  is  transmitted 
P(  1/0)  =  probability  of  a  “1”  being  received  when  a  “0”  is  transmitted 

For  this  research  effort,  a  binary  phase  shift  keying  (BPSK)  modulation  scheme 
will  be  assumed  due  to  the  coherent  optical  detection.  In  addition,  antipodal  signals 
will  be  assumed.  The  P(e)  for  antipodal  BPSK  is  written  (34:158) 


where 


Q(-)  =  complementary  error  function  =  f£°  exp  (— y)  du 

Eb  =  signal  energy  per  binary  symbol  (J) 

N0  =  additive  white  Gaussian  noise  (AWGN)  energy  (J) 


(60) 


For  coherent  optical  systems,  the  noise  is  dominated  by  shot  noise  (which 
has  Poisson  statistics)  rather  than  thermal  noise  (which  has  Gaussian  statistics) 
(28:2158);  nevertheless,  a  Gaussian  distribution  can  be  assumed  for  the  noise  (33:431). 
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The  quantity  Eb/N0  is  related  to  the  SNR  by  (34:158) 


E±^S_B , 
N0  Nt  R 

where 


(61) 


S/Nr  =  total  SNR 
S  =  average  signal  power  ( W  •  m2) 
Nt  =  average  noise  power  (W  •  m2) 
Ba  =  signal  bandwidth  (Hz) 

R  =  bit  rate  (bps) 


The  total  SNR,  S/Nt ,  must  now  be  related  to  the  normalized  SNR  derived  in 
Section  3.3.2  (Equation  (58)).  First,  the  average  noise  power,  Nt ,  is  related  to  N 
(noise  per  unit  bandwidth)  by 

Nt  =  NxBs  (62) 


Second,  the  relationship  between  the  average  signal  power,  S',  and  the  average  signal 
power  per  watt  of  cw  output  power,  can  be  expressed  as 


S  = 


<s> 

WA 


Pi 


(63) 


where  Pt  is  the  output  power  of  the  transmitter  laser.  The  transmitter  (laser)  power, 
Pt,  is  related  to  the  effective  cw  output  power  by 


Pt  = 


Wa 

D\ 


(64) 


Examining  the  expression  for  <  S  >,  Equation  (26),  from  a  dimensional  anal¬ 
ysis  viewpoint,  it  can  be  seen  that  <  S  >  is  expressed  in  units  of  W  •  m4.  Noting 
that  Nt  is  expressed  in  units  of  W  -  m2,  then  ^7  is  expressed  in  units  of  m2.  Thus, 
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this  expression  for  the  average  SNR  must  be  normalized  by  to-2.  Recall  that  W a, 
the  normalization  factor  used  in  developing  the  normalized  SNR  for  this  study,  is 
defined  as  the  intensity  of  the  transmitter  field  integrated  over  the  transmitter  aper¬ 
ture  (that  is,  Wa  is  a  power  density)  (28:2158).  Dividing  the  average  SNR  by  the 
aperture  dimension  will  then  restore  the  dimensionless  integrity  of  the  SNR.  That 
is, 

S_  =  <£>  Wd  f65^ 

Nt  NtWa  D\  K  1 

Equations  (62)  and  (65)  can  be  combined  to  yield 


S  _<S>  WA 
Nt  ~  NWa  BsD2t 


so  that  substituting  Equations  (61)  and  (66)  into  Equation  (60)  results  in  the  BER 
equation  to  be  evaluated  for  this  feasibility  study: 


P(e)  =  Q 


<S>  WA  \* 
NWa  RD't  ) 


(67) 


3.4  Summary 

In  this  chapter,  the  general  SNR  was  developed.  Under  the  assumptions  of 
uniformly  weighted  apertures  and  a  Gaussian  cross-correlation  function  for  the  tur¬ 
bulence,  an  expression  for  a  closed  form  solution  to  the  general  SNR  was  obtained. 
Also,  under  the  assumption  of  square  apertures,  a  different  solution  to  the  SNR  was 
obtained,  this  one  being  solved  using  numerical  integration  methods.  Finally,  an 
expression  for  the  BER  was  derived  based  on  the  BERs  relationship  to  the  SNR. 
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IV.  Analysis  and  Evaluation 


This  chapter  addresses  the  analysis  and  evaluation  of  the  SNR  and  BER  data 
gathered  during  this  research  effort.  Specifically,  this  chapter  addresses  five  topics: 
trends  analysis  of  SNR  versus  aperture  sizes,  SNR  efficiency  factor,  SNR  versus 
altitude  and  aperture  sizes,  C„  comparison,  and  BER  feasibility  analysis.  The  trends 
analysis  of  the  SNR  versus  aperture  dimensions  investigates  the  relationship  between 
the  SNR  and  the  effective  transmitter  and  receiver  diameters  when  assuming  infinite 
Gaussian  weighting  functions.  In  addition,  the  relationship  between  the  SNR  and 
the  phase  coherence  length,  p0,  is  investigated.  The  SNR  efficiency  factor  analysis 
discusses  the  SNR  degradation  due  to  turbulence.  The  SNR  versus  altitude  and 
aperture  sizes  analysis  compares  the  SNRs  for  various  altitudes  as  well  as  the  SNRs 
for  various  aperture  sizes.  The  comparison  compares  the  SNRs  determined  with 
the  Hufnagel  model  (23)  with  Feldmann’s  (10)  experimentally  derived  C\  data.  The 
BER  feasibility  analysis  discusses  the  BER  data  and  presents  the  evaluation  of  the 
feasibility  analysis. 

4-1  Signal-to-Noise  Ratio  Versus  Aperture  Size  Trends  Analysis 

In  this  section,  the  results  of  Lutomirski  and  Buser’s  closed  form  solution  for  the 
SNR  is  investigated.  Recall  that  the  closed  form  solution  was  obtained  by  assuming 
infinite,  circularly  symmetric,  Gaussian  transmitter  and  receiver  aperture  weighting 
functions. 

The  reason  for  investigating  the  infinite  Gaussian  weighting  functions  SNR  is 
two  fold.  First,  the  general  behavior  of  the  SNR  with  respect  to  varying  the  diameter 
of  the  apertures  is  evaluated.  This  behavior  will  assist  in  determining  the  expected 
behavior  of  the  square  aperture  SNR.  Second,  the  effect  of  wavelength  on  the  SNR 
is  analyzed. 
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4.1.1  Expected  Results  The  behavior  of  the  infinite  Gaussian  weighting 
functions  SNR  was  investigated  by  Lutomirski  and  Buser.  In  their  investigation, 
they  found  that  the  SNR  improves  with  increasing  transmitter  and  receiver  aperture 
diameters  until  each  diameter  is  on  the  order  of  the  phase  coherence  length,  p0 - 
Once  the  aperture  diameters  exceeded  p0)  the  SNR  degraded  so  that  any  further 
increase  in  the  aperture  diameters  does  not  result  in  a  corresponding  increase  in  the 
SNR  (28:2159).  In  other  words,  the  SNR  will  saturate  once  the  aperture  diameters 
increase  beyond  p0.  (Recall  Fried’s  and  Chen’s  saturation  discussion  from  Section 
2.3.4.) 

The  other  result  expected  from  this  section  concerns  the  behavior  of  the  SNR 
with  respect  to  the  wavelength,  A.  Lutomirski  and  Buser  pointed  out  that,  while 
the  SNR  for  propagation  in  a  vacuum  is  independent  of  wavelength,  the  turbulence 
affected  SNR  is  not.  They  discussed  the  fact  that,  since  p0  increases  with  wavelength, 
potentially  better  performance  can  be  expected  for  longer  wavelengths  (28:2159). 

4-1.2  Varying  Aperture  Sizes  Results  Table  7  lists  the  parameters  used 
for  this  particular  analysis.  In  addition,  the  phase  coherence  length,  p0>  is  provided 
for  comparative  purposes.  Figures  10  through  12  illustrate  the  results  of  Equa¬ 
tion  (39).  (Note  that,  for  these  figures,  the  physical  receiver  aperture  diameter  is 
assumed  to  equal  the  effective  receiver  aperture  diameter,  and  W  represents  W^.) 

Figures  10  through  12  show  the  normalized  SNR,  in  dB,  with  respect  to  the 
effective  aperture  diameter,  d.  Figure  10  illustrates  the  case  where  the  SNR  has 
reached  saturation  for  all  aperture  diameters  since,  from  Table  7,  p0  is  on  the  order 
of  micrometers.  Figures  11  and  12  show  similar  results  for  the  z  =  100  km  and 
z  =  160  km  data  points  (in  which  p0  is  slightly  less  than  1  cm).  However,  for  the 
z  —  40  km  data  points,  p0  is  equal  to,  or  slightly  greater  than,  1  cm.  In  this  case, 
it  can  be  seen  that  the  SNR  saturates  at  approximately  3  dB  above  the  SNR  for 
d  =  1  cm.  This  is  in  agreement  with  Fried’s  and  Chen’s  findings  as  discussed  in 
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(ap)  MN/<S> 


Table  7.  Infinite  Gaussian  Aperture  Signal-to-Noise  Ratio  Parameters 


Figure 

A  (nm) 

tSSBBU 

z(km) 

p0{cm) 

40 

10 

1.5 

1.625E-15 

100 

160 

2.31E-4 

40 

MEM 

11 

830 

1.625E-15 

mm 

o 

|g| 

■9 

40 

1.1 

12 

904 

1.625E-15 

100 

0.7 

160 

0.5 

Figure  10.  Infinite  Gaussian  Aperture  Signal-to-Noise  Ratio:  A  =  1.5  nm 
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<S>/NW  (dB)  ^  <S>/NW  (dB) 


’’igure  11.  Infinite  Gaussian  Aperture  Signal-to-Noise  Ratio:  A  =  830  nm 


Figure  12.  Infinite  Gaussian  Aperture  Signal-to-Noise  Ratio:  A  =  904  nm 
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Section  2.3.4.  Thus,  with  respect  to  varying  aperture  sizes,  the  results  obtained  are 
as  expected. 

4-1.3  Wavelength  Evaluation  Results  Referring  to  Figures  10  through 
12,  the  results  with  respect  to  wavelength  are  also  as  expected.  Specifically,  Fig¬ 
ures  10  and  11  depict  the  significantly  better  performance  of  A  =  830  nm  and 
A  =  904  nm  than  the  A  =  1.5  nm  performance  depicted  by  Figure  10.  As  a  conse¬ 
quence  of  these  results,  and  noting  that  the  performance  of  the  904  nm  wavelength 
is  slightly  better  than  830  nm,  only  the  904  nm  wavelength  will  be  considered  for 
the  remainder  of  this  study. 

4.1.4  Summary  The  behavior  of  the  infinite  Gaussian  weighting  functions 
SNR,  as  derived  by  Lutomirski  and  Buser,  was  investigated.  It  was  shown  that 
the  derived  SNR  performed  as  expected  with  respect  to  increasing  aperture  sizes, 
saturation  behavior  and  wavelength. 

4.2  Signal-to-Noise  Ratio  Efficiency  Factor  Analysis 

The  remainder  of  this  chapter  addresses  the  SNR  expression  that  was  derived 
with  the  assumptions  of  uniformly  weighted,  finite,  square  apertures.  In  this  section, 
a  SNR  performance  measure  is  analyzed.  Specifically,  a  SNR  efficiency  factor,  a,  is 
defined  to  be  the  ratio  of  the  SNR  with  turbulence  effects  considered  to  the  SNR 
without  turbulence  effects  considered.  This  factor,  a,  is  thus  one  measure  of  the  SNR 
degradation  due  to  turbulence.  The  SNR  efficiency  analysis  is  performed  to  verify 
that  the  derived  SNR  expression  using  uniformly  weighted,  finite,  square  apertures 
behaves  in  a  manner  consistent  with  Fried’s  (infinite  plane  wave  incident  on  a  finite 
receiver  aperture)  and  Lutomirski  and  Buser’s  (infinite  Gaussian  transmitter  and 
receiver  aperture  weighting  functions)  results. 

In  particular,  this  section  will  briefly  address  the  overall  efficiency  of  the  uni¬ 
formly  weighted,  finite,  square  apertures  SNR.  This  overall  efficiency  is  interesting 
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Table  8.  Best,  Medium  and  Worst  Case  Signal-to-Noise  Ratio  Efficiencies 
Parameters 


A  (nm) 

D(cm) 

V(m/s) 

h(m ) 

best 

904 

1 

218 

5000 

medium 

904 

7.5 

218 

7500 

worst 

904 

15 

218 

10000 

primarily  in  that  it  provides  a  baseline  for  the  remaining  discussions.  The  behavior 
of  a  with  respect  to  varying  aperture  sizes  and  altitudes  is  then  examined.  Finally, 
the  relationship  between  a  and  the  normalized  SNR  is  discussed. 

4.2.1  Expected  Results  As  briefly  discussed  in  Section  4.1.1,  the  SNR 
performance  is  expected  to  degrade  when  a  circular  aperture  diameter  exceeds  the 
phase  coherence  length,  p0.  This  behavior  can  easily  be  extended  to  the  case  of 
square  apertures.  That  is,  the  SNR  performance  is  expected  to  degrade  when  the 
square  aperture  side  length  exceeds  p0.  Also,  consistent  with  Fried’s  and  Chen’s 
saturation  discussion,  it  is  expected  that,  as  a  decreases,  the  SNR  approaches  the 
saturation  limit  (for  the  particular  parameters  used). 

4.2.2  Overall  Efficiency  Figure  13  shows  the  best,  medium  and  worst  ef¬ 
ficiencies  calculated  for  a  wavelength  of  904  nm.  As  can  be  seen,  by  varying  the 
aperture  dimension,  D ,  and  the  altitude,  h,  the  efficiency,  a,  ranges  from  0.16  to 

67.2  percent.  Table  8  contains  the  parameters  used  in  plotting  Figure  13. 

The  importance  of  Figure  13  is  best  understood  when  compared  to  the  appli¬ 
cable  phase  coherence  lengths.  Table  9  shows  the  value  of  p0  for  various  altitudes 
and  path  lengths.  By  comparing  these  values  with  Figure  13,  it  becomes  apparent 
why  the  best  case  (D  =  1  cm)  results  in  such  a  significantly  better  a  than  the  other 
two  cases  (D  =  7.5  cm,  D  =  15  cm).  The  phase  coherence  lengths  are  less  than  1 
cm  for  the  medium  and  worst  case  efficiencies;  whereas,  the  phase  coherence  length 
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Figure  13.  Best,  Medium  and  Worst  Case  SNR  Efficiencies:  A  =  904  nm 


for  the  best  case  is  3.40  cm. 

Thus,  the  behavior  of  a  is  as  predicted  in  Section  4.2.1.  It  is  seen  that,  the 
more  the  square  aperture  side  length  exceeds  p0)  the  worst  the  efficiency,  a. 


4.2.3  Alpha  Versus  Varying  Parameters  This  section  examines  the  be¬ 
havior  of  a  while  varying  different  parameters  as  detailed  in  Table  10.  The  results 


Table  9.  Phase  Coherence  Lengths:  A  =  904  nm,  V  =  218  m/s 


h(m) 

5000 

7500 

10000 

12500 

z[km) 

Po(cm) 

Po{cm) 

p0{cm) 

p0{cm) 

40 

3.40 

1.37 

1.09 

1.28 

60 

2.67 

1.07 

0.86 

1.01 

80 

2.24 

0.90 

0.72 

0.85 

100 

1.96 

0.79 

0.63 

0.74 

120 

1.76 

0.71 

0.56 

0.66 

140 

1.60 

0.64 

0.51 

0.60 

160 

1.48 

0.60 

0.48 

0.56 

60 


Table  10.  Alpha  Versus  Varying  Parameters  Cases 


Figure 

D(cm) 

h(m) 

z(-km) 

mm 

varying 

Bfl 

varying 

16 

varying 

12500 

mm 

varying 

5000-12500 

40 

B9 

varying 

5000-12500 

100 

K9 

5000-12500 

160 

20 

varying 

21 

|D| 

varying 

22 

15 

.  varying 

are  presented  graphically  in  Figures  14  through  22.  First,  the  behavior  of  a  with 
respect  to  aperture  size  is  examined.  Figures  14  through  19  illustrate  this  behav¬ 
ior  for  both  increasing  path  length  (Figures  14  through  16)  and  increasing  altitude 
(Figures  17  through  19). 

As  can  be  seen  from  Figures  14  through  16,  and  recalling  the  pQ  values  defined 
in  Table  9,  a  decreases  with  increasing  aperture  side  length,  as  expected.  However, 
it  is  interesting  to  note  that  the  relative  decrease  in  a  is  less  at  the  longer  path 
lengths,  z.  This  is  due  to  the  dependence  of  p0  on  z  (37:1403): 

p0  =  (o.545fc2C^z)~*  (68) 


That  is,  as  z  increases,  p0  decreases. 

Figures  17  through  19  also  illustrate  the  behavior  of  a  with  increasing  aperture 
side  length.  What  is  interesting  to  note  here,  however,  is  the  behavior  of  a  with 
respect  to  increasing  altitude,  h  (for  any  given  aperture  side  length).  It  is  seen 
that  a  does  not  continually  decrease  as  h  increases,  but  rather  reaches  a  minimum 
at  about  10000  m  and  then  starts  increasing  again.  This  behavior,  however,  is  as 
expected  when  considering  the  behavior  of  the  turbulence  strength,  C2.  For  all  of 
these  figures,  C2  was  modeled  using  the  Hufnagel  model.  As  illustrated  in  Figure  23, 
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Figure  14.  Alpha  Versus  Path  Length  with  Varying  Aperture  Siz e:h  =  5000  m 


Figure  15.  Alpha  Versus  Path  Length  with  Varying  Aperture  Size:/i  =  8750  m 
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Figure  16.  Alpha  Versus  Path  Length  with  Varying  Aperture  Siz e:h  =  12500  m 


Figure  17.  Alpha  Versus  Altitude  with  Varying  Aperture  Size:z  =  40  km 
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alpha  alpha 


Figure  19.  Alpha  Versus  Altitude  with  Varying  Aperture  Size:z  =  160  km 
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the  Hufnagel  model  characterizes  the  turbulence  strength  as  increasing  over  5000  m 
to  about  10000  m,  then  decreasing  to  12500  ra.  (A  good  discussion  on  the  behavior 
of  C *  can  be  found  in  the  Infrared  Handbook  (23),  Chapter  6.) 

The  behavior  of  a  with  respect  to  varying  altitude  is  further  illustrated  in 
Figures  20  through  22.  Again,  this  behavior  is  as  expected. 

4.2.4  Alpha  Versus  Turbulence  Affected  Signal-to-Noise  Ratio  As  a 
final  evaluation  of  the  behavior  of  a,  this  section  examines  the  relationship  between 
a  and  the  normalized,  turbulence  affected  SNR.  To  reiterate,  the  purpose  for  this 
evaluation  is  to  verify  that  the  uniformly  weighted,  square  aperture  SNR  derived  in 
this  study  produces  results  consistent  with  previous  studies.  Table  11  presents  the 
values  of  a  and  the  normalized  SNR  for  varying  aperture  dimensions  and  altitudes, 
while  maintaining  a  constant  path  length.  Specifically,  Table  11  investigates  the  case 
for  A  =  904  nm,  V  —  218  m/s,  and  z  =  100  km. 

As  can  be  seen,  Table  11  illustrates  the  expected  behavior  of  both  a  and  the 
normalized  SNR  with  respect  to  p0.  That  is,  as  the  aperture  side  length,  D,  is 
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Figure  22.  Alpha  Versus  Path  Length  with  Varying  Altitude:  D  =  15  cm 
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Figure  23.  Plot  of  Turbulence  Strength  Characteristics  Using  the  Hufnagel  Model 


Table  11.  Alpha  Versus  Signal- to-Noise  Ratio:  A  =  904  nm, 
V  =  218  m/s ,  z  —  100  km 


h(m)  p0(cm) 


5000  1.96 


increased  beyond  the  phase  coherence  length,  p0,  a  decreases  while  the  normalized 
SNR  begins  to  saturate.  Note  that,  if  the  normalized  SNR  data  were  extrapolated 
for  the  case  of  D  =  p0,  it  appears  that  the  saturated  SNR  is  about  3  dB  greater  (as 
predicted  by  Chen  (5)). 

4-2.5  Summary  The  behavior  of  the  SNR  efficiency  factor,  a,  was  exam¬ 
ined.  It  was  shown  that  a  behaved  as  expected  in  all  respects.  It  was  further  shown 
that  the  behavior  of  a  versus  the  normalized  SNR  supports  the  findings  of  Fried  (11) 
and  Chen  (5)  with  respect  to  a  SNR  saturation  limit. 

4.3  Signal-to-Noise  Ratio  Versus  Altitude  and  Aperture  Sizes 

In  the  previous  section,  a  versus  varying  altitudes  and  aperture  sizes  was  exam¬ 
ined.  In  this  section,  attention  is  focused  on  the  behavior  of  the  turbulence  affected 
SNR  while  varying  the  altitude  and  aperture  sizes.  This  investigation  of  the  turbu¬ 
lence  affected  SNR  is  accomplished  for  the  purpose  of  again  verifying  the  consistency 
of  this  study’s  results  with  previous  results  (such  as  Fried’s  (11)  and  Lutomirski  and 
Buser’s  (28)). 

4.3. 1  Expected  Results  The  results  expected  from  the  uniformly  weighted, 
square  aperture  SNR  versus  altitude  and  aperture  size  are  the  same  as  the  expected 
results  of  Section  4.2.1.  To  reiterate,  as  the  aperture  side  length,  D ,  is  increased 
beyond  the  phase  coherence  length,  p0,  the  SNR  is  expected  to  reach  a  saturation 
limit  (which  is  about  3  dB  greater  than  the  SNR  achieved  for  D  =  p0).  In  addition, 
as  the  altitude  is  increased,  the  SNR  is  expected  to  behave  in  a  manner  inversely 
related  to  the  turbulence  strength,  C%. 

4.3.2  Signal-to-Noise  Ratio  Versus  Altitude  First,  the  behavior  of  the 
SNR  with  respect  to  altitude  is  examined.  Table  12  lists  the  parameters  used  for 
Figures  24  through  27.  These  figures  illustrate  the  behavior  of  the  turbulence  affected 
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Table  12.  Signal-to-Noise  Ratio  Versus  Altitude  Parameters 


Figure 

A(nm) 

V(m/s) 

D(cm) 

z(km) 

WA{mW) 

24 

904 

218 

1 

Varying 

1 

25 

904 

218 

Varying 

1 

26 

904 

218 

Varying 

1 

27 

904 

218 

Varying 

100 

1 

h  (m) 


Figure  24.  Signal-to-Noise  Ratio  Versus  Altitude  with  Varying  Path  Length  - 
D  =  1  cm 


SNR  with  respect  to  altitude.  As  discussed  in  Section  4.2.2,  the  behavior  shown  in 
Figures  24  through  27  is  as  expected.  That  is,  as  the  altitude  increases  from  5000  m 
to  10000  m,  the  SNR  decreases;  whereas,  as  the  altitude  increases  further  to  12500  m, 
the  SNR  increases.  Note  that  this  is  true  regardless  of  path  length  and  ai  <•  -  ,-ire  size. 


4-3.3  Signal-to-Noise  Ratio  Versus  Aperture  Size  The  behavior  of  the 
SNR  with  respect  to  aperture  size  is  examined  here.  Figures  26  through  27  also 
illustrate  the  behavior  of  the  turbulence  affected  SNR  with  respect  to  aperture  size. 
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<S>/N  (dB)  5  <S>/N  (dB) 


25.  Signal-to-Noise  Ratio  Versus  Altitude  with  Varying  Path  Length  - 
D  =  10  cm 


Figure  26.  Signal-to-Noise  Ratio  Versus  Altitude  with  Varying  Aperture  Size  - 
z  =  40  km 
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Figure  27.  Signal-to-Noise  Ratio  Versus  Altitude  with  Varying  Aperture  Size  - 
z  =  100  km 

As  can  be  seen,  these  figures  support  the  findings  of  Section  4.2.2.  It  is  also  inter¬ 
esting  to  observe  the  behavior  of  the  SNR  with  respect  to  aperture  side  length  in 
a  format  similar  to  Fried’s  and  Chen’s  saturation  limit  figure,  Figure  3.  Figure  28 
shows  the  normalized  SNR  versus  D/ p0  for  A  =  904  nm,  h  =  7500  m,  z  =  100  fcm, 
and  p0  =  0.79  cm.  As  can  be  seen,  Figure  28  presents  a  very  close  approximation  to 
Figure  3. 

4.3.4  Summary  The  behavior  of  the  turbulence  affected  SNR  for  uniformly 
weighted,  finite,  square  apertures  was  examined.  It  was  shown  that  the  SNR  per¬ 
formed  as  expected  for  both  cases  considered:  varying  altitude  and  varying  aperture 
sizes. 


4.4  Signal-to-Noise  Ratio  Versus  C\  Comparison 

In  this  section,  the  SNRs  determined  using  the  Hufnagel  model  (23)  for  C\  are 
compared  to  the  SNRs  determined  using  Feldmann’s  experimentally  derived  C\  (10). 
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Table  13.  C\  Comparison  (10:55,58) 


Constants 

D  =  10.01cm 

WA  =  l.Om  W 

A  =  904 nm 

ms^rn 

Run 

HESuflHi 

mmssmm 

v\m/t)  | 

Feldmann 

Hufnagel 

Feldmann 

Hufnagel 

1 

9754 

5.0B-13 

1.2B-15 

114.55 

145.86 

2 

9144 

198.1 

5.0E-13 

1.4E)-15 

114.55 

3 

2.5E-13 

4.1E-17 

118.15 

159.54 

4 

3353 

154.3 

5.0E-14 

4.1E-17 

135.62 

166.28 

5 

6096 

46300 

5.0E-14 

3.8E-16 

126.53 

151.54 

6 

12497 

46300 

159.5  | 

5.0E-14 

7.2B-16 

126.53 

148.41 

This  comparison  is  made  in  an  effort  to  determine  the  impact  on  the  SNR  when  a 
theoretical  model  is  used  to  determine  C%  as  opposed  to  experimentally  derived  C\. 
Table  13  lists  the  parameters  used  and  presents  the  resulting  SNRs.  As  can  be  sej 
Feldmann’s  values  are  100  to  1000  times  greater  than  the  predicted  by  . 
Hufnagel  model;  that  is,  Feldmann  predicts  much  stronger  turbulence  for  all  altitude?; 
than  the  Hufnagel  model.  As  a  result,  the  SNRs  determined  with  Feldmann’s  C \ 
are  much  smaller. 

Several  factors  should  be  considered  when  evaluating  the  C\  comparison  re¬ 
sults.  First,  the  Hufnagel  model  is  just  that,  a  model.  It  does  not  accurately  match 
the  observed  fine  structure  of  the  turbulence,  but  it  is  more  accurate  than  other 
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models  used  to  predict  C*  (23:6-15).  Second,  C\  is  a  continuously  varying  feature 
of  turbulence.  For  example,  wind  sheer,  time  of  day  and  temperature  inversions  all 
affect  the  turbulence  strength  (23:6-13  -  6-14).  Finally,  Feldmann’s  data  is  based 
on  a  “best-fit”  approach.  As  stated  by  Feldmann,  “...values  of  C\  were  consistently 
larger  than  predicted  by  the  Hufnagel  model  but,  for  the  most  part  remained  within 
the  accepted  range  of  variability”  (10:62). 

4.5  Bit  Error  Rate  Analysis 

In  this  section,  the  heart  of  this  study  is  discussed:  feasibility.  First,  the  BER 
data  collected  is  analyzed  with  respect  to  effective  cw  output  power  and  data  rate. 
Then,  the  feasibility  of  an  air-to-air  optical  heterodyne  communication  system  is 
evaluated.  However,  before  this  analysis  is  presented,  several  concerns  need  to  be 
addressed.  First,  it  should  be  emphasized  again  that  this  analysis  only  considers  the 
losses  due  to  propagation  path  length  and  the  turbulence-induced  phase  perturba¬ 
tions.  As  will  be  discussed  in  more  detail  shortly,  the  losses  due  to  beam  steering, 
beam  spreading,  and  scintillation  are  ignored.  Also,  the  transmitter  and  receiver, 
as  well  as  receiver  tracking,  are  assumed  to  be  ideal.  Further,  as  was  shown  in  Sec¬ 
tion  4.4,  the  Hufnagel  model  for  C\  predicts  SNRs  consistently  larger  (on  the  order 
of  30  dB)  than  Feldmann’s  experimentally  derived  data.  All  of  these  concerns 
serve  to  reduce  the  required  effective  cw  catput  power,  and  thus  the  transmitter 
(laser)  power,  necessary  to  achieve  the  desired  BERs.  This  will  be  addressed  further 
in  Section  4.5.3. 

4.5.1  Expected  Results  Consistent  with  the  SNR  results,  the  BERs  are 
expected  to  “saturate”  for  aperture  side  lengths,  D ,  greater  than  p0.  That  is,  the 
ranges  at  which  some  defined  BER  is  achieved  should  not  increase  appreciably  as  D 
increases  for  D  >  p0. 
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4-5.2  Bit  Error  Rate  Data  Analysis  The  BER  data  will  follow  the  same 
trends  as  the  turbulence  affected  SNR.  This  is  obvious  upon  examination  of  the 
relationship  between  the  BER  and  SNR: 


BER  =  Q 


\SNR 


l  R  J 


(69) 


where 


<?(•)  =  complementary  error  function 
R  =  data  rate  (bps) 

It  is  also  obvious  from  Equation  (69)  that  as  the  data  rate,  R,  increases,  the  BER 
decreases. 

Tables  14  through  17  present  examples  of  these  trends.  In  these  tables,  the 
data  points  presented  represent  the  range  (km)  at  which  the  appropriate  data  rate 
is  obtained. 

4-5.3  Air-to-Air  Optical  Heterodyne  Communication  System  Feasi¬ 
bility  The  results  of  this  analysis  demonstrate  that  an  air-to-air  optical  heterodyne 
communication  system  is  feasible.  However,  this  study  did  not  consider  many  of  the 
additional  losses  that  would  affect  the  communication  link  and  the  system,  such  as 
beam  steering,  beam  spreading,  scintillation  and  the  receiver  and  transmitter  perfor¬ 
mance  characteristics.  All  of  these  considerations  will,  of  course,  reduce  the  overall 
SNR  and  consequently  the  BER.  Further  reductions  in  the  SNR  and  BER  may  be 
realized  when  more  accurate  turbulence  strength,  C%,  factors  are  considered.  The 
question  becomes,  when  the  additional  loss  considerations  are  taken  into  account, 
will  acceptable  BERs  still  be  obtained  for  the  ranges  presented  in  this  study. 

To  answer  this  question,  the  effective  cw  transmitter  output  powers,  Wa,  used 
to  generate  Tables  14  through  17  need  to  be  examined.  Recall  that  Wa  is  related  to 
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Table  14.  Bit  Error  Rate  Trends  -  R  =  20000  bps:  -  Range  greater  than 

160  km ;  -  Range  less  than  40  km 


BER 

—3  ^ 

- w=n 

h(m ) 

5000 

7500 

10000 

12500 

5000 

7500 

10000 

12500 

D(cm) 

WA{nW) 

Pt{nW) 

0.1 

mm 

+ 

+ 

mm 

+ 

+ 

+ 

152 

+ 

i 

0.5 

EH 

+ 

+ 

Bg§ 

+ 

+ 

+ 

+ 

+ 

1.0 

10000 

+ 

+ 

KB 

+ 

+ 

+ 

+ 

+ 

5.0 

50000 

+ 

+ 

KH 

+ 

+ 

+ 

+ 

+ 

■■i 

0.1 

17.78 

+ 

+  ' 

+ 

+ 

+ 

+ 

+  • 

+ 

0.5 

88.89 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

BB 

1.0 

177.78 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

1  1 

5.0 

888.89 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

0.1 

4.44 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

15 

0.5 

22.22 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

1.0 

44.44 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

5.0 

222.22 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

Table  15.  Bit  Error  Rate  Trends  -  R  =  40000  bps:  -  Range  greater  than 

160  km ;  -  Range  less  than  40  km 


BER 

-—3 

- IF*5 

-■ 

h(m) 

5000 

7500 

10000 

12500 

5000 

7500 

10000 

12500 

D{cm) 

WA(nW) 

Pt(nW) 

- 

0.1 

+ 

+ 

158 

+ 

+ 

130 

118 

127 

1 

0.5 

l 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

1.0 

1 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

J 

5.0 

| 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

“ 

1 

■Sf 

0.1 

17.78 

+ 

+ 

+ 

+ 

+ 

155 

135 

149 

0.5 

88.89 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

}\ 

BH 

1.0 

177.78 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

i 

i 

^B 

5.0 

888.89 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

■■ 

0.1 

4.44 

+ 

4- 

i 

+ 

BHI 

+ 

155 

135 

149 

15 

0.5 

22.22 

+ 

+ 

+ 

Bfl 

+ 

+ 

+ 

+ 

/j 

fi 

1.0 

44.44 

+ 

+ 

-  + 

mm 

+ 

+ 

+ 

+ 

5.0 

222.22 

+ 

+ 

+ 

KB 

+ 

+ 

+ 

+ 

4 
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Table  16.  Bit  Error  Rate  Trends  -  R  =  100000  bps:  -  Range  greater  than 

160  km\  -  Range  less  than  40  km 


BER 

IF5  “ 

- 

h(m) 

5000 

7500 

10000 

12500 

5000 

7500 

10000 

12500 

D(cm) 

WA(nW) 

0.1 

158 

125 

T— 1 

i-H 

122 

109 

91 

84 

89 

i 

0.5 

w  • --'in 

+ 

+ 

+ 

+ 

+ 

+ 

152 

+ 

1.0 

10000 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

5.0 

50000 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

psi 

0.1 

17.78 

+ 

149- 

130 

144 

+ 

116 

101  • 

112 

0.5 

88.89 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

I 

1.0 

177.78 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

5.0 

888.89 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

0.1 

4.44 

+ 

150 

130 

144 

+ 

117 

101 

112 

15 

0.5 

22.22 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

1.0 

44.44 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

5.0 

222.22 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

Table  17.  Bit  Error  Rate  Trends  -  R  =  1000000  bps:  -  Range  greater  than  j 

160  km]  -  Range  less  than  40  km  } 


BER 

nr3 

lO-io 

h(m ) 

5000 

7500 

10000 

12500 

5000 

7500 

10000 

12500 

D(cm) 

WA{nW) 

Pt(nW ) 

0.1 

53 

48 

46 

47 

- 

- 

- 

- 

1 

0.5 

.  l 

114 

95 

88 

93 

78 

68 

64 

67 

1.0 

l 

158 

125 

114 

122 

109 

91 

84 

89 

5.0 

+ 

+ 

+ 

+ 

+ 

+ 

152 

+ 

■ 

0.1 

17.78 

126 

72 

63 

70 

97 

56 

49 

54 

0.5 

88.89 

+ 

120 

104 

115 

+ 

93 

81 

89 

■ 

1.0 

177.78 

+ 

149 

130 

144 

+ 

116 

101 

112 

i 

5.0 

888.89 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

0.1 

4.44 

128 

73 

63 

70 

100 

56 

54 

15 

0.5 

22.22 

+ 

121 

105 

116 

+ 

94 

90 

1.0 

44.44 

+ 

150 

130 

144 

+ 

117 

101 

112 

5.0 

222.22 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

76 


Table  18.  Transmitter  Powers 


Pt  (nW) 

WA  (nVF) 

0.1 

0.5 

1.0 

5.0 

H 

1000 

5000 

10000 

50000 

D  (cm) 

17.78 

88.89 

177.78 

888.89 

15 

4.44 

22.22 

44.44 

222.22 

the  transmitter  power,  Pt ,  by 

WA  =  p,d2t 

(where  transmitter  power  refers  to  the  intensity  of  the  transmitter  field  incident  on 
the  transmitting  aperture).  Table  18  consolidates  the  associated  transmitter  powers 
for  the  effective  cw  transmitter  output  powers  used  to  generate  Tables  14  through 
17. 

The  transmitter  powers  presented  here  are  much  lower  than  one  would  expect 
for  the  application  at  hand.  Transmitter  (laser)  powers  on  the  order  of  hundreds  of 
milliwatts  would  seem  more  reasonable.  But,  again,  this  analysis  was  accomplished 
while  only  considering  the  losses  due  to  propagation  path  length  and  turbulence 
induced  phase  perturbations.  The  SNR  efficiency  analysis  in  Section  4.2  showed  that 
the  turbulence  induced  phase  perturbations  caused  the  SNR  efficiency  to  be  reduced 
as  much  as  99.9  percent.  If  the  losses  due  to  beam  steering,  beam  spreading,  and 
scintillation  also  produced  significant  reductions  to  the  SNR  efficiency,  plus  taking 
into  account  less  than  ideal  transmitter  and  receiver  performance  (including  tracking 
by  the  receiver)  as  well  as  stronger  turbulence  (that  is,  greater  C*),  it  becomes 
evident  that  the  overall  efficiency  of  the  optical  communication  system  may  decrease 
significantly.  Thus,  the  transmitter  power  would  need  to  be  increased  considerably 
to  account  for  the  additional  losses.  For  example,  assume  the  additional  losses 
combine  to  reduce  the  efficiency  by  another  99.9  percent  and  the  combined  non¬ 
ideal  transmitter  and  receiver  performance  (including  tracking)  contributes  another 
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99.9  percent  efficiency  loss.  The  transmitter  (laser)  power  would  then  need  to  be 
increased  on  the  order  of  104  to  get  similar  results  to  those  presented  in  Tables  14 
through  17. 

4-6  Summary 

This  chapter  presented  the  analysis  and  evaluation  of  the  SNR  and  BER  data 
gathered  during  this  study.  It  was  shown  that  the  SNR  efficiency  ranged  from  0.16 
to  67.2  percent  for  the  parameters  considered  in  this  study.  The  behavior  of  the 
SNR  efficiency  factor,  a,  was  examined  with  respect  to  altitude,  path  length  and 
aperture  size.  As  expected,  a  decreases  with  increasing  path  length  and  increasing 
aperture  size  ( D  >  p0),  and  is  inversely  related  to  the  turbulence  strength, 

The  turbulence  affected  SNR  was  then  investigated.  It  also  exhibited  the  expected 
behavior  for  increasing  path  length  (SNR  decreased),  aperture  size  (SNR  saturated 
for  D  >  p0 )  and  turbulence  strength  (SNR  inversely  related).  A  comparison  of  the 
C\  model  used  for  this  study  was  compared  to  experimentally  derived  C\  data  from 
a  SNR  standpoint.  From  this  comparison  it  was  shown  that  the  model  predicts 
weaker  turbulence  than  the  experimentally  derived  data,  thus  the  SNRs  were  higher 
for  the  model.  Finally,  the  BER  performance  was  presented  with  respect  to  data  rate 
and  path  length  (ranges).  This  analysis  illustrated  that  BERs  up  to  10-10  could  be 
achieved  at  ranges  of  160  km  under  certain  conditions.  The  final  result  of  the  BER 
analysis  was  that  air-to-air  optical  heterodyne  communication  systems  are  feasible. 
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V.  Conclusions  and  Recommendations 

This  chapter  presents  the  conclusions  of  this  research  effort.  Specifically,  the 
findings  of  the  BER  feasibility  study  are  presented.  Following  the  conclusions,  rec¬ 
ommendations  for  further  research  axe  discussed. 

5.1  Conclusions 

The  foremost  question  to  be  answered  by  this  thesis  effort  was  “is  an  air-to-air 
optical  heterodyne  communcation  system  feasible?”.  From  a  BER  standpoint,  the 
answer  is  a  resounding  “yes”!  The  analysis  in  Section  4.5  illustrates  that  (depending 
on  the  system  parameters)  communications  can  be  conducted  at  ranges  in  excess  of 
160  km  for  BERs  of  10~10  and  data  rates  of  20  to  1000  kbps.  Of  course,  even  longer 
ranges  are  possible  for  a  BER  of  10-5 . 

In  addition,  it  was  shown  that  aperture  dimensions,  D,  greater  than  the  phase 
coherence  length,  p0,  result  in  a  saturation  of  the  SNR.  This  saturation  is  approx¬ 
imately  4.5  dB  greater  than  the  SNR  achieved  when  D  =  p0.  Therefore,  the  best 
choice  for  the  aperture  dimensions  are  on  the  order  of  the  phase  coherence  length. 

5.2  Recommendations 

Given  that  air-to-air  optical  communication  systems  are  feasible,  the  ques¬ 
tion  now  becomes,  as  discussed  in  Section  4.5,  what  transmitter  (laser)  powers  are 
required  to  achieve  these  results.  As  discussed  thro  ghout  this  study,  this  analy¬ 
sis  assumed  nearly  ideal  conditions.  That  is,  this  study  only  considered  the  losses 
due  to  path  length  propagation  and  turbulence  induced  phase  perturbations.  This 
leads  to  several  possible  research  topics  in  this  area  of  air-to-air  optical  heterodyne 
communication  systems. 
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1.  Investigate  the  effects  of  additional  losses  on  the  system,  to  include  beam 
steering,  beam  spreading,  scintillation,  transmitter  characteristics  (such 
as  laser  phase  noise)  and  receiver  characteristics  (such  as  quantum  effi¬ 
ciency,  shot  noise  and  tracking  capability). 

2.  Compare  this  system  to  a  direct  detection  system.  As  discussed  in  Chap¬ 
ter  II,  the  coherent  detection  system  performs  better  when  turbulence 
is  not  a  factor  (that  is,  in  a  vacuum).  How  much  more  the  turbulence 
degrades  the  coherent  detection  system  than  the  direct  detection  system, 
and  the  overall  impact  on  the  absolute  performance  would  better  enable 
the  Air  Force  in  determining  which  system  is  more  feasible  from  a  cost 
versus  performance  standpoint. 

3.  Experientally  verify  the  results  of  the  research,  or  at  least  some  portion  of 
the  results.  This,  of  course,  is  the  ideal  method  to  proving  feasibility,  but 
also  the  most  costly.  It  may  be  possible  to  verify  portions  of  the  research 
which,  in  turn,  would  lead  to  better  theoretical  results. 

4.  Investigate  the  use  of  adaptive  optics  schemes  to  compensate  for  the  wave 
front  tilt  induced  by  the  turbulence.  Included  in  this  effort  may  be  an 
investigation  into  tracking  techniques  by  the  receiver. 

5.  Investigate  the  use  of  an  array  of  small  receivers,  each  on  the  order  of 
the  phase  coherence  length,  as  a  means  of  increasing  the  efficiency  of  the 
communication  system. 
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Appendix  A.  Phase  Coherence  Length  Discussion 

In  this  appendix,  various  representations  for  the  phase  cohere^,.  length  are 
briefly  discussed.  Specifically,  five  different  nomenclature  and  mathematical  defini¬ 
tions  used  to  represent  the  phase  coherence  length  are  presented.  It  is  then  demon¬ 
strated  that  each  of  these  representations  produce  approximately  the  same  values. 

A.l  Phase  Coherence  Length  Nomenclature  and  Definitions 

In  this  section,  the  nomenclatures  and  mathematical  definitions  of  five  separate 
papers  are  presented:  Fried  (11),  Lutomirski  and  Buser  (28),  Yura  (37),  Greenwood 
and  Fried  (21),  and  Wilkins  (36). 

Fried  refers  to  the  phase  coherence  length  as  the  “efficiency  saturation  dimen¬ 
sion”  and  denotes  it  r0.  He  stated  that,  for  horizontal  propagation,  r0  can  be  defined 
as  (11:57,63) 

r0  =  (0.1089A-2C2z)-3/5  (70) 

where 

k  =  optical  wave  number  (equal  to  ^)  (to-1) 
z  =  propagation  path  length  (m) 

C\  —  turbulence  strength  (to-2/3) 

Lutomirski  and  Buser  refer  to  the  phase  coherence  length  as  both  a  “turbulence- 
induced  coherence  length”  (28:2153)  and  a  “transverse  coherence  length”  (28:2159). 
They  denote  the  phase  coherence  length  by  p0  and  define  it  to  be  (28:2157) 

p0  =  (0.5 k2C*z)  3/5  (71) 
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Yura  refers  to  the  phase  coherence  length,  p0 ,  as  “...the  MCF  of  a  point  source 
located  in  the  plane  of  observation  and  observed  in  the  aperture  plane”  (37:1403). 
He  then  defines  p0  for  homogeneous  turbulence  by  (37:1403) 

p0  =  (0.545fc2C2z)“3/5  (72) 


Greenwood  and  Fried  refer  to  the  phase  coherence  length  as  “Fried’s  coherence 
length”  and  denote  it  r0.  They  define  r0  by  (21:200) 

t0  =  (o.423&2  C\z )  "3/S  s-1  (73) 


where 


for  plane  waves 
for  spherical  waves 


and  where  z'  is  defined  to  be  an  observation  point  along  the  propagation  path  such 
that  0  <  f  <  1. 


Finally,  Wilkins  refers  to  the  phase  coherence  length  as  the  “atmospheric  phase 
coherence  length”  and  denotes  it  r0.  He  also  stated  that  r0  is  known  as  “the  diffrac¬ 
tion  limited  aperture  of  the  atmosphere”  or  as  “Fried’s  parameter”  (36:1).  Wilkins 
defines  r0  by  the  expression  (36:1) 


0.423&2  l*  C*(L)dL 
.  Jo 


-3/S 


(74) 


where  L  defines  the  propagation  path  integral  (m).  By  modeling  C2  using  the 
Hufnagel  model,  C2  is  not  a  function  of  the  path  length  z.  Thus,  Equation  (74)  can 
be  rewritten  as 

t0  =  (o.423fc2C2z)  ^  (75) 

which  is  equivalent  to  Equation  (73)  for  s  =  1. 
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Table  19.  Phase  Coherence  Length  Comparison  Parameters 


Example 

A  (nm) 

Clim- -*/*) 

z{km\ 

1 

1.5 

1.353  x  10"lS 

40 

2 

830 

1.625  X  10“15 

100 

3 

904 

2.661  x  10"16 

160 

Table  20.  Phase  Coherence  Length  Comparison 


Example 

Fried 
r0(cm ) 

Lutomirski /Buser 
p0(cm) 

Yura 

p0{cm) 

Greenwood/Fried  ' 
r0(cm ) 

Wilkins 

r0(cm) 

1 

7.09  x  10"4 

6.22  X  10"4 

5.91  x  10"4 

6.88  x  10"4 

6.88  x  10"4 

2 

0.72 

0.63 

0.60 

0.70 

0.70 

3 

1.77 

,1.58 

1.48 

1.72 

1.72 

A.2  Phase  Coherence  Length  Comparisons 

In  this  section,  the  different  phase  coherence  length  representations  are  com¬ 
pared.  It  is  seen  that  these  representations  all  predict  approximately  the  same  values 
for  the  phase  coherence  length. 

Table  19  presents  the  parameters  used  to  determine  the  phase  coherence  lengths 
as  defined  by  Equations  (70)  through  (75).  The  phase  coherence  lengths  are  pre¬ 
sented  in  Table  20. 
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Appendix  B.  Signal-to-Noise  Ratio  for  Circular  Apertures 


In  this  appendix,  the  SNR  for  uniformly  weighted,  circular  apertures  will  be 
developed.  Starting  with  the  general  SNR  equation  developed  in  Chapter  III,  a 
Gaussian  cross- correlation  function  will  be  assumed.  Sum  and  difference  variable 
substitutions  will  then  be  made,  followed  by  reduction  of  the  integral  equation  to  a 
four-fold  expression. 

The  uniformly  weighted,  circular  aperture  SNR,  denoted  SNRC ,  is  derived 
from  the  general  SNR,  Equation  (30).  Restated  here  for  convenience,  the  general 
SNR  is  written 


<S> 

nwa 


£(;)  sSf?  JJJJnwswto*'  W 

xH  (2,  ■-rs)<P&<Pra<iJ&«lJ& 


(76) 


where 

V'A  (£,)  n  fe)  =  ^§-Ua  {u)U'a  fo)  (77) 

and  where  the  terms  are  defined  in  Chapter  III. 

For  uniformly  weighted,  circular  apertures,  the  complex  disturbance  at  r, 
Ua  (l)>  can  be  written 

UA(r)  =  U0W(r)  (78) 

and  the  receiver  weighting  factor,  F ^(g),  can  be  written 

F(t)  =  W  (1)  (79) 
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where 


W(x)  = 


i  iii<f 

o  fel  >  7 


and  where  D  =  aperture  diameter. 

Recall  that  Wa,  the  normalization  factor  is  defined  to  be 


wA  =  J 

Substituting  Equations  (78)  and  (80)  into  (81)  then  yields 

WA  =  J  J  !/;H'(rl)r(E!)4,fe 

r  rDr/2  „ 

=  J  J-DT/,U°dL'dL* 

=  VlD\ 


(80) 


(81) 


(82) 


where  Dt  =  transmitter  aperture  diameter. 

The  uniformly  weighted,  circular  aperture  SNR,  SNRC ,  can  now  be  written  by 
substituting  Equations  (77),  (78),  (79)  and  (82)  into  (76): 


SNRC  = 


xH  (&  _  2.V &  _  &)  <^hd2r2d2lld2l2 


(83) 


To  be  consistent  with  the  development  of  the  uniformly  weighted,  square  aper¬ 
tures  SNR  development  of  Chapter  III,  a  Gaussian  cross-correlation  function  is  again 
assumed.  Thus,  SNRC  can  now  be  written 
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xexp{i 

X  exp  {“7  [l(&  “&)“(&“  E«)|2]}  <84) 

Note  that  |  (p^  —  g2)  —  (rj  —  r2)|2  can  also  be  written  |  (j^  —  —  (g2  -  r2)  |  . 

In  order  to  reduce  Equation  (84),  sum  and  difference  variable  substitutions  are 
made  as  follows: 

Ax,-  =  £.  -  r 

=  \  (li  +  Li)  (85) 

so  that 

p.  =  Sx,-  +  xAx.- 

t,  — «  2 — , 

U  =  Sa-iAa  (86) 

Substituting  Equations  (85)  and  (86)  into  (84),  and  rearranging  the  integrals 
then  yields 

SNR-  -  2^p//eXP|l  ^ ‘|!  ”  l4$2|2l} 

X  exp  (— i  [|Aii  -  Afel2]  | 
x  {//  W  (&£,  -  (&r2  -  jAjj) 

Ax?^  d2Sx1d2 Sx? 

xd2Ax1d2Ax?  (87) 


x  W  (exx  +  ^Axj)  W  (Sx2  +  i 
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Figure  29.  Circular  Aperture  Geometry  (21:205) 

Looking  at  the  inner  integrals: 

Area  =  J  W  (y,xx  -  ^Agi)  W  (jfa  +  ^Ax,) 

x  J  w  (sx2  -  W  (sx2  +  ^Ax^  d2 Sx2  (88) 

This  area  is  illustrated  in  Figure  29.  The  limits  of  integration  are  described  by  the 
area  of  overlap  of  the  circles.  By  defining  a  new  variable  I  to  be  (21:196) 

I<(1)  =  Jw(pxi-  ^Ax,-)  W  (p x{  +  iAx,-)  d2£x,.  (89) 


where 


Ar  =  area  of  receiver  aperture 
At  =  area  of  transmitter  aperture 


the  inner  integrals  can  now  be  written  by  substituting  Equation  (89)  into  (88)  to 
yield 

Area  =  A^I,(1)I2(1)  (90) 


Greenwood  solves  the  integral  J<,  Equation  (89),  with  the  result  (21:196) 


m = 


j  cos-1  (Axi/Dx)  —  (Axi/Dj)  1  —  (Ax;/jDx)2j 


1/2 


0  <  Axi  <  Dj 
Df  <  A Xi 


(91) 


Thus,  Equation  (88)  can  now  be  rewritten 

Area  =  |cos-1  (Axi/Dt)  -  (Axi/Dt)  [l  —  (Aa^/Dy)2]  ^  j 

x  {cos"1  (A x2/Dt)  -  (A*j/Pr)  jl  -  (A x2fDT)2} 1/2 }  (92) 


The  final  form  of  the  uniformly  weighted,  circular  aperture  SNR  can  now  be 
written  by  substituting  Equation  (92)  into  (87)  and  by  making  use  of  the  relation 
Ar  =  ttjD^/4: 


x  jcos-1  (Axi/Dt)  —  {Axi/Dt)  [i  ~  {Axi/Dt)2]  ^  j 
x  |cos_1  (AX2/ Dt)  —  (A^/jDt)  [l  —  (A12 / jDj*)2]  ^  j 
X^As^As?  (93) 

Note  that  Equation  (93)  can  not  be  separated  due  to  the  cross  term  introduced  by 

exp  {-^  [|A*i  -  Ax2|2]}- 
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Appendix  C.  Signal-to-Noise  Ratio  Computer  Code 

Development 

The  signal-to-noise  ratio  (SNR)  developed  in  Chapter  III  (Equation  (58))  must 
be  solved  using  numerical  methods.  This  appendix  details  the  development  of  the 
numerical  methods  used,  the  manipulation  of  the  SNR  equation  into  a  codable  format 
and  a  listing  of  the  FORTRAN  77  computer  code. 

C.l  Numerical  Integration  Methods 

The  development  of  the  numerical  method  to  be  used  is  based  on  the  approach 
of  “...adding  up  the  value  of  the  integrand  of  a  sequence  of  abscissas  within  the 
range  of  integration”  (31:103).  Figure  30  illustrates  the  geometry  of  the  approach. 
Mathematically,  the  geometry  of  Figure  30  can,  be  expressed  as 

x;  =  x0  +  ih  (94) 

/(* f)  =  U  (95) 


where 


i  —  0, 1,2, ...,  N  +  1 
/(•)  —  integrand 

h  =  constant  step  size  within  the  abscissa  sequence 


A  method  to  solve  such  a  sequence  is  the  trapezoidal  rule.  The  trapezoidal 
rule  can  be  written  as  (31:105) 


jT  /(*)<&= &  1 + !/>] + 0  (k3r) 


(96) 


where 
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Figure  30.  Numerical  Integration  Geometry.  Numerical  methods  formulas  with 
equally  spaced  abscissas  compute  the  integral  of  a  function  between 
x0  and  ajjv+i.  Closed  formulas  evaluate  the  function  on  the  boundary 
points,  while  open  formulas  refrain  from  doing  so  (useful  if  the  evaluation 
algorithm  breaks  down  on  the  boundary  points).  (31:104) 

f"  —  second  derivative  of  the  function  f(x) 

O(-)  =  error  term 


The  error  term,  0(»),  is  undetermined,  otherwise  the  integral  would  be  solvable  in  a 
closed  form. 

Equation  (96)  is  a  two-point  formula.  Press  and  others  show  that  this  can  be 
extended  to  a  three-point  formula  by  use~of  Simpson’s  rule.  Simpson’s  rule  can  be 
written  (31:105) 


[Xi 

/  f(x)dx  =  h 
Jx  1 


llf' +  5/!  +  lh 


+  o(as/(4)) 


(97) 


Press  and  others  further  show  that,  by  using  Equation  (96)  N  —  1  times  (that 
is,  for  each  of  the  intervals)  and  adding  up  the  results,  an  extended  trapezoidal  rule 
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can  be  obtained.  The  extended  trapezoidal  rule  is  written  (31:107) 


h  2^  +  /a  +  /a+ 

- f  fN- 1  + 


(98) 


By  applying  the  extended  trapezoidal  rule,  Equation  (98),  to  Simpson’s  rule, 
Equation  (97),  the  extended  Simpson’s  rule  results.  The  extended  Simpson’s  rule 
can  be  written  (31:108) 


txs 

/  f{x)dx  =  h 
Jx  l 


liA  +  5/2  +  + 1/<+ 

*  +  3/^-2  +  -/n-1  +  jj/fV 


+ 


°Gft) 


(99) 


Press  and  others  use  the  extended  Simpson’s  rule,  along  with  the  extended 
trapezoidal  rule,  to  write  subroutines  to  calculate  an  integral  numerically;  These 
subroutines  (QSIMP  and  TRAPZD)  are  used  to  numerically  solve  the  inner  integrals 
of  the  two-fold  integrals  of  the  SNR  equation.  The  outer  integrals  are  then  solved 
numerically  by  again  applying  the  extended  Simpson’s  rule. 

Once  the  SNR  is  calculated,  the  BER  is  determined.  As  discussed  in  Chapter 
III,  the  BER  uses  the  complimentary  error  function,  £?(•).  Another  form  of  the 
complimentary  error  function,  erfc,  is  written 

(100) 


erfc(x)  =  *4=  J  cxp  (-ti2)  du 


where  (34:743) 

Q{x)  =  \erfc(j*) 


Press  and  others  write  an  erfc  routine,  which  in  turn  uses  other  routines,  which  are 
also  employed  in  this  code.  The  routines  are  ERFC,  GAMMP,  GAMMQ,  GAMMLN, 
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GSER,  and  GCF  (31:157-164). 


C.2  Signal-to-Noise  Ratio  Equation  Manipulation  into  a  Codable  Format 

The  SNR  equation  to  be  modeled,  Equation  (58),  presents  itself  in  an  easily 
codable  format  with  two  exceptions.  The  first  “problem”  concerns  the  imaginary 
term  in  the  exponential: 

ik 


exp 


2z 


(AxJj  -  A®’3) 


where  q  denotes  either  ®  or  y.  The  second-  “problem”  arises  due  to  the  nature  of  the 
three  functions  being  integrated: 


D  -  |A®,,-| 


exp 


exp 


l—  (a*;,  -  a*j2) 
1 


L  P‘ 


(Ax,i  -  A®,2) 


A  solution  to  the  imaginary  term  “problem”  is  straight  forward.  First,  separate 
the  variable  terms,  then  apply  Euler’s  identities.  Thus 


exp 


Yz  (A*?1  -  A*y 


=  exp 


exp 


ik  .  2 


and 


exp 


4^ 


In  this  form,  it  can  easily  be  verified  that  the  imaginary  terms  cancel  out  since 
the  integration  limits  on  Ax,i  and  A®,2  are  the  same.  The.result  of  this  manipulation 
is 


exp 


Yz  (Aa&  “  Aa&) 


=>•  cos 


(Aa^)  cos  ^A*;2)+sin 


,\  .  tk 
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The  second  problem  requires  an  understanding  of  the  behavior  of  the  three 
terms  being  integrated.  Figure  31  illustrates  the  case  for  Axqi  =  0.10  while  varying 
Axq2  over  the  range  of  integration.  For  the  figure,  A  =  904  nm,  z  —  100  km  and 
p0  =  0.6  cm.  Note  that 


tri 

■■4 

D  - 

cos /sin 

—4 

cos  | 

H 

— y 

exp 

i 

( 


k  N 

—At2 

2z'glj 


cos 


'  u  > 
—Ax2 
,2  *  ^ 


+  sin 


(  AaxJ.)  sin  ( 


—Ax2  N 
2  z*X'\ 


-—(Axql-Axq2) 


As  can  be  seen  from  Figure  31,  the  H  term  significantly  limits  the  area  of  integration. 
To  compensate  for  this,  the  limits  of  integration  for  the  model  are  set  to  5cr  where 
a  is  defined  as  the  variance  for  the  Gaussian  function: 
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C.S  FORTRAN  77  Code 


The  following  is  the  listing  of  the  computer  code  used  to  determine  the  SNR 
along  with  the  associated  BER.  The  code  is  written  in  the  FORTRAN  77  computer 
language. 
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C  PROGRAM:  snr.f 

C  AUTHOR:  Capt  Rebecca  N.  Seeger,  GE-90D 

C  This  program  calculates  the  Signal-to-Noise  Ratio  for 
C  an  optical  heterodyne  communication  system  in  the 
C  presence  of  atmospheric  turbulence.  The  only  losses 

C  considered  are  those  associated  vith  the  path  length 

C  and  those  associated  vith  the  phase  perturbations 

C  resulting  from  propagation  through  an  inhomogeneous 

C  medium  (i.e.  turbulence). 

C  This  code  calculates  the  SNR  assuming  both  transmitter 

C  and  receiver  are  at  the  same  altitude  (>  3000  m) .  The 

C  SNR  is  run  for  path  lengths  of  40,  45,  ...  160  km. 

C  Variable  Definitions: 

C  STEP  -  the  number  of  segments  that  the  outer  integral 

C  is  broken  up  into  when  applying  Simpson's  rule. 

C  DIST  -  the  number  of  path  lengths  considered. 

C  IAREA#  -  inner  integral  function  to  be  integrated 

C  using  the  routine  QSIMP. 

C  S#  -  the  value  of  the  integration  returned  from  QSIMP. 

C  0AREA#  -  outer  integral  function  to  be  integrated 

C  using  the  routine  INTEGRATE. 

C  SO#  -  the  value  of  the  integration  returned  from 

C  INTEGRATE. 

C  X#  -  the  integration  variables  (representing  transmitter 

C  and  receiver  square  aperture  dimensions) . 

C  LIM#  -  the  limits  of  integration  for  the  turbulence 

C  affected  SNR. 

C  Z,Z0  -  path  length. 

C  RH0,RH00  -  phase  coherence  length. 

C  C  -  constant  in  front  of  the  integrals. 

C  CN2  -  turbulence  strength  as  defined  by  the  Hufnagel 

C  model . 

C  CN21,CN22  -  the  two  terms  in  the  Hufnagel  model. 

C  D1  -  transmitter  aperture  dimension. 

C  D2  -  receiver  aperture  dimension. 

C  LAMBDA  -  wavelength. 

C  ETA  -  quantum  efficiency  of  the  receiver. 

C  K  -  optical  wave  number  (equal  to  2*PI/LAMBDA) . 

C  E  -  charge  on  an  electron  (equal  to  1.602e-19). 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PI  -  3.14159265359 

V  -  velocity  of  the  aircraft  (relative  wind  speed) . 

H  -  altitude  of  the  aircraft. 

SIGMA  -  standard  deviation  of  the  cross-correlation 
term  in  the  integration.  This  defines  the 
limits  of  integration  for  the  turbulence 
affected  SNR. 

W  -  effective  cw  output  power  (or  power  density) . 

R  -  data  rate. 

SNRTRB  -  turbulence  affected  SNR  for  purpose  of 

determining  ALPHA.  Does  not  contain  terms 
in  common  with  SNRCLN. 

SNRCLN  -  non-turbulence  affected  SNR  (see  above 
with  respect  to  SNRTRB) . 

ERROR  -  argument  of  the  erfc  routine  used  to  determine 
the  BER. 

BER  -  bit  error  rate. 

ALPHA  -  ratio  of  SNRTRB  to  SNRCLN. 

SNR  -  actual  turbulence  affected  SNR.  Output 

value  is  expressed  as  a  ratio  (not  in  dB) . 

I#,J#  -  do  loop  counters. 

NINT,NREAL  -  number  of  integration  steps  used  in  the 
INTEGRATE  routine. 

RESULTS  -  File  results  will  be  written  to. 


C  ************************************************************ 


C  Variable  declarations. 


INTEGER  STEP,  DIST 
PARAMETER  (STEP=401) 
PARAMETER  (DIST=25) 


C  Integration  variables. 

REAL  IAREA1,  IAREA2,  SI,  S2,  IAREA3,  IAREA4,  S3,  S4 
REAL  0AREA1 ( 1 : STEP) ,  0AREA2 ( 1 : STEP) ,  OAREA ( 1 : STEP) 
REAL  0AREA3 ( 1 : STEP) ,  0AREA4 ( 1 : STEP) 

REAL  SOI,  S02,  S03,  S04,  XI,  X2,  LIM1,  LIM2 


C  Constant  variables. 
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REAL  Z(1:DIST) ,  RH0(1:DIST),  C 

REAL  CN2,  CN2i,  CN22,  Dl,  D2,  LAMBDA,  ETA 

REAL  ZO,  RHOO,  K,  E,  PI,  V,  H,  SIGMA,  W,  R 


C  Output  variables . 

REAL  SNRTRB(1:DIST) ,  BER(1:DIST),  ERROR 
REAL  SNRCLN(1:DIST) ,  ALPHA(1:DIST) ,SNR(i:DIST) 

C  Other  variables. 

INTEGER  II,  12,  13,  Jl,  J2,  33,  NINT 
INTEGER  J4,  J5 
REAL  NREAL 
CHARACTER+20  RESULTS 

C  Declare  data,  functions. 

EXTERNAL  I AREA 1,  IAREA2,  ERFC,  GAMMP,  GAMMQ,  GAMMLN 
EXTERNAL  IAREA3,  IAREA4 

INTRINSIC  COS,  SIN,  EXP,  ABS,  SQRT,  MIN,  MAX 
DATA  PI/3.14159265359/,  E/1.602E-19/ 

C  Declare  common  variables. 

COMMON/BLOCKl/Xl , ZO , RHOO , K , Dl , D2 , SIGMA 
C0MM0N/BL0CK2/0AREA 

C  Input  user  defined  terms. 


PRINT*, ’Enter  results  filename:’ 
READ  (*, ’ (A20) ’)  RESULTS 


PRINT*, 'Enter  transmitter  aperture  side  length  in  cm:’ 
READ*,D1 


PRINT*, ’Enter  receiver  aperture  side  length  in  cm:’ 
READ*,D2 

PRINT*, ’Enter  wavelength  in  nm:’ 

READ*, LAMBDA 

PRINT*, ’Enter  quantum  efficiency:’ 
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HEAD*, ETA 


PRINT* , ’  Enter  aircraft  altitude  in  m:’ 

READ*,H 

PRINT* , ’Enter  relative  wind  speed  in  m/s:’ 
READ*,V 

PRINT*, ’Enter  transmitter  output  power  in  nW:’ 
READ*,W 

PRINT*,  ’Enter  data  rate  in  bps:’ 

READ*,  R 


C  Calculate  the  remaining  constants. 

Di  =  D1  *  '1 .OE-2 
D2  =  D2  *  1 . OE-2 

CN21  =  8.2E-30  *  (V**2)  *  (H**5)  *  (H**5)  *  EXP(-H/1000.0) 
CN22  =  2.7E-16  *  EXP(-H/1500.0) 

CN2  =  (CN21  *  1.0E-26)  +  CN22 

LAMBDA  =  LAMBDA  *  1.0E-9 
K  =  (2.0  *  PI) /LAMBDA 

C  =  ( (K**2) *ETA) / (8 . 0* (PI**2) * (D2**2) *E) 

W  =  W  *  1.0E-9 

C  The  SNR  is  calculated  over  the  path 

C  lengths  between  40  and  160  km. 

DO  10  II  =  1,DIST 

Z(I1)  =  (5.0  *  REAL(Il-l)  +40.0)  *  1.0E3 
RHO(Ii)  =  (0.545  *  (K**2)  *  CN2  *  Z(Il))**(-0.6) 

10  CONTINUE 

C  Determine  the  SNR. 

DO  20  12=1, DIST 
ZO  =  Z(I2) 
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RHOO  =  RH0(I2) 

SIGMA  =  (5.0  *  RH00)/SQRT(2.0) 

DO  30  13*1, STEP 

XI  *  -Dl+(2.0  *  D1  *  (REAL (13-1) /REAL (STEP- 1))) 
IF((ABS(D1-ABS(X1)) .LT.1.0E-6) .OR. (ABS(D2-ABS(X2)) . 

+  LT.1.0E-6))  THEN 

51  *  0.0 

52  =  0.0 
ELSE 

LIM1  *  MAX(X1-SIGMA,X1-D2) 

LIM2  *  MIN (Xl+SIGMA ,Xi+D2) 

C  Calculate  the  inner  integrals  over  XI. 

CALL  QSIMP(IAREA1,LIM1,LIM2,S1) 

CALL  QSIMP (IAREA2 , LIM1 , LIM2 , S2) 

CALL  QSIMP (IAREA3,-D2,D2, S3) 

CALL  QSIMP (IAREA4,-D2,D2,S4) 

END  IF 

0AREA1(I3)  *  SI* (Dl-ABS (Xl) ) *COS ( (K* (Xl**2) ) / (2 . 0*Z0) ) 
0AREA2(I3)  *  S2* (Dl-ABS (XI) ) *SIN ( (K* (Xi**2) ) / (2 . 0*Z0) ) 
0AREA3(I3)  *  S3* (Dl-ABS (XI) ) *COS ( (K* (Xl**2) ) / (2 . 0*Z0) ) 
0AREA4(I3)  *  S4* (Dl-ABS (XI) )*SIN ( (K* (Xl**2) ) / (2 . 0*Z0) ) 
30  CONTINUE 

C  Now  calculate  the  outer  integrals  over  X2. 

NREAL  =  REAL(STEP-l)/2.0 
NINT  *  INT (NREAL) 

DO  40  J 1=1, STEP 

OAREA(Jl)  =  OAREAl(Jl) 

40  CONTINUE 

CALL  INTEGRATE(S01 ,NINT,-D1 ,D1) 

DO  45  J2=i,STEP 

0AREA(J2)  =  0AREA2(J2) 

45  CONTINUE 

CALL  INTEGRATE (S02, NINT, -Di ,D1) 

DO  47  J4=l .STEP 
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0AREA(J4)  =  0AREA3(J4) 

47  CONTINUE 

CALL  INTEGRATE (S03,NINT,-D1 ,D1) 

DO  49  J5*i,STEP 

0AREA(J5)  *  0AREA4CJ5) 

49  CONTINUE 

CALL  INTEGRATE(S04,NINT,-D1,D1) 

C  Determine  the  SNR  values. 

SNRTRB(I2)  *  (SOI  +  S02)**2 
SNRCLN(I2)  =  (S03  +  S04)**2 
ALPHA (12)  =  SNRTRB (12)/ SNRCLN (12) 

SNR(I2)  *  (W*C*SNRTRB (12) ) / ( (Z (12) **2) * (Dl**2) * (D2**2) ) 

20  CONTINUE 

C  Determine  the  BER  using  the  erfc  function,  then  convert 

C  to  Q  function. 

DO  50  J3=1,DIST 

ERROR  =  SQRT(SNR(J3)/R) 

BER(J3)  =  0.5  *  ERFC (ERROR) 

50  CONTINUE 

C  Print  out  the  results  to  a  user  named  file. 

OPEN (UNIT=20 ,FILE=RESULTS , STATUS* »NEH » , 

+  ACCESS*' SEQUENTIAL' ,F0RM*' FORMATTED' , 

+  ERR*1000) 

DO  72  J3=l,7 

WRITE (20 , 73 , ERR=2000)  '  ' 

73  FORMAT(A) 

72  CONTINUE 

WRITE(20,76,ERR*2000) 'Di*' ,D1, 'D2=' ,D2, 'LAMBDA*' .LAMBDA, 

+  'H*',H,'C*',C 

76  F0RMAT(7X,A,F5.3,4X,A,F5.3,4X,A,EU.4,4X,A,F8.0,4X,A,E11.4) 


WRITE(20,88,ERR=2000) 'W=' ,W, 'R=' ,R, ' LTA* ' ,LiA,'V=' ,V> 'Cn2=' ,CN2 


100 


88 


F0RMAT(7X,A,E11.4,4X,A,F8.0,4X,A,F6.3,4X,A,F8.0,4X,A,E11.4) 


WRITE(20,73,ERR«2000)  *  .  * 

WRITE(20,73,ERR«2000)  »  • 

DO  60  J3*1,DIST 

WRITE (20 ,75, ERR*2000)  ' Z« ' , Z ( J3) , ' SNR (v/) • ' , SNRTRB ( J3) , 

+  'SNR(w/o)»' ,SNRCLN(J3) , 'ALPHA-' ,ALPHA(J3) 

+  ,'SNR«',SNR(J3),'BER-',BER(J3) 

75  F0RMAT(7X,A,F8.0,4X,A,E11.4,4X,A,E11.4,4X,A,E11.4,4X,A, 

+  Ell. 4, 4X, A, Eli. 4) 

60  CONTINUE 

ENDFILE  20 

CLOSE  (UNIT*20 ,  STATUS88 » KEEP ' ) 

GO  TO  999 

1000  PRINT*, ’CANNOT  OPEN  FILE! * ,ERRSIG 

GO  TO  999 

2000  PRINT*, 'ERROR  WRITING  TO  FILE! ' ,ERRSIG 

GO  TO  999 

999  END 

C  Define  the  inner  integral  functions. 

REAL  FUNCTION  IAREA1(X2) 

COMMON/BLOCKl/Xi ,Z0 ,RHOO ,K,D1 ,D2 , SIGMA 
REAL  X1,X2,Z0,RH00,D2,K 

I AREA 1  =  (D2-ABS(X2»  *  C03((K*(X2**2))/(2.0*Z0))  * 

+  EXP ( (- 1 . 0/RH00**2) * ( (X1-X2) **2) ) 

RETURN 

END 


REAL  FUNCTION  IAREA2(X2) 

REAL  Xi,X2,D2,K,Z0,RH00 
COMMON/BLOCKl/Xi , ZO , RHOO ,K , D1 ,D2 , SIGMA 
IAREA2  =  (D2-ABS(X2))  *  SIN((K*(X2**2))/(2.0*Z0))  * 
+  EXP( (-1 . 0/RH00**2)*( (X1-X2) **2) ) 

RETURN 
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REAL  FUNCTION  IAREA3(X2) 

REAL  X2,D2,K,Z0 

COMMON/BLOCKl/Xi ,ZO ,RH00 ,K,D1 ,D2 .SIGMA 
IAREA3* (D2-ABS (X2)  ) *COS ( (K* (X2**2) ) / (2 . 0*Z0) ) 
RETURN 
END 

REAL  FUNCTION  IAREA4(X2) 

REAL  X2.D2.K.Z0 

COMMON/BLOCKl/Xi ,ZO,RHOO ,K,Di ,D2 .SIGMA 
IAREA4* (D2-ABS (X2)  ) *SIN ( (K* (X2**2) ) / (2 . 0*Z0) ) 
RETURN 
END 


Outer  integral  routines. 


SUBROUTINE  INTEGRATE (AREA , N , A , B) 

C0MM0N/BL0CK2/0AREA 

REAL  AREA,  A,  B,  HI,  OAREA(lOOl) 

INTEGER  N,  IFACTOR,  L 
HI  =  (B-A)/ (2.0  *  REAL(N)) 

IFACTOR  =  4 
AREA  =0.0 
DO  50  L=2,2*N 

AREA  =  AREA  +  REAL (IFACTOR)  *  OAREA(L) 
IFACTOR  =  8/IFACTOR 
CONTINUE 

AREA  =  AREA  +  OAREA(i) 

AREA  =  AREA  +  0AREA(2*N+i) 

AREA  =  AREA  *  (Hl/3.0) 

RETURN 

END 

Inner  integral  routines. 


SUBROUTINE  QSIMP(FUNC , A*B ,S) 
PARAMETER  (EPS=5.E-3,  JMAX=20) 
0ST=-1 . E30 


0S«  -1.E30 
DO  11  J**l,  JMAX 

CALL  TRAPZD(FUNC,A,B,ST, J) 

S* (4 . *ST-0ST)/3 . 

IF  (ABS(S-OS).LT.EPS*ABS(OS))  RETURN 
OS*S 
0ST*ST 
11  CONTINUE 

PAUSE  ’Too  many  steps.' 

END 


SUBROUTINE  TRAPZD(FUNC,A,B,S,N) 
IF  (N.EQ.l)  THEN 
S=0 . 5* (B-A) * (FUNC (A) +FUNC (B) ) 
IT«1 
ELSE 
TNM=IT 

DEL=(B-A)/TNM 
X*A+0 . 5*DEL 
SUMsO . 

DO  11  J*1,IT 
SUM*SUM+FUNC(X) 

X=X+DEL 
11  CONTINUE 

S=0. 5* (S+ (B-A) *SUM/TNM) 
IT=2*IT 
ENDIF 
RETURN 
END 

C  BER  routines. 

FUNCTION  ERFC(X) 

IF(X.LT.O. )THEN 
ERFC=1 . +GAMMP ( . 5 , X**2) 

ELSE 

ERFC=GAMMQ(.5,X**2) 

ENDIF 

RETURN 

END 
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FUNCTION  GAMMP(A.X) 

IF (X . LT . 0 . .OR.A.LE.O.)PAUSE 
IF(X.LT. A+1.)THEN 
CALL  GSER (GAMSER ,  A ,  X , GLN) 

GAMMP«GAMSER 

ELSE 

CALL  GCFCGAMMCF,A,X,GLN) 

GAMMP-1  ..^GAMMCF 
ENDIF 
RETURN 
END 

FUNCTION  GAMMQ(A,X) 

IFCX.LT. 0. .OR.A.LE.O. )PAUSE 
IF(X.LT.A+1.)THEN 
CALL  GSERCGAMSER,A,X,GLN) 

GAMMQ-1. -GAMSER 
ELSE 

CALL  GCFCGAMMCF.A.X.GLN) 

GAHMQ'GAHMCF 

ENDIF 

RETURN 

END 

FUNCTION  GAMMLNCXX) 

REAL*8  C0FC6) »STP, HALF, ONE, FPF,X,TMP,SER 

DATA  COF , STP/76 . 18009173D0 , -86 . 50532033D0 , 24 . 01409822D0 , 

*  - 1 . 231739516D0 , . 120858003D-2 , - . 536382D-5 , 2 . 50662827465D0/ 

DATA  HALF , ONE, FPF/O . 5D0 , 1 . ODO , 5 . 5D0/ 

X=XX-ONE 

TMP*X+FPF 

TOP*  CX+HALF)*LOG  CTMP)  -TMP 
SER*ONE 
DO  11  >1,6 
X*X+ONE 

SER=SER+COF  C J) /X 
11  CONTINUE 

GAMMLN*TOP+LOG CSTP*SER) 

RETURN 

END 

SUBROUTINE  GCF  CGAMMCF , A , X , GLN) 
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PARAMETER  (ITMAX-100, EPS-3. E-7) 
GLN-GAMMLN(A) 

GOLD-O. 

AO-1. 

Al-X 

BO-O. 

Bl-1. 

FAC-1. 

DO  11  N«l, ITMAX 
AN-FLOAT(N) 

ANA-AN-A 

AO- (Al+AO+ANA) *FAC 
BO- (Bl+BO+ANA) *FAC 
ANF-AN+FAC 
Al»X*AO+ANF*Al 
B1-X*B0+ANF*B1 
IF(Al.NE.O.)THEN 
FAC-1. /A1 
G-B1+FAC 

IF(ABS((G-GOLD)/G) .LT.EPS)CD  TO  1 
GOLD-G 
ENDIF 

11  CONTINUE 

PAUSE  'A  too  large,  ITMAX  too  small' 

1  GAMMCF-EXP (-X+A*ALOG (X) -GLN) *G 

RETURN 
END 

SUBROUTINE  GSER(GAMSER,A,X,GLN) 
PARAMETER  (ITMAX-lOO, EPS-3. E-7) 
GLN-GAMMLN (A) 

IF(X.LE.O. )THEN 
IF (X.LT.O.) PAUSE 
GAMSER-0 . 

RETURN 
ENDIF 
AP-A 
SUM-1. /A 
DEL-SUM 

DO  11  N-l, ITMAX 
AP-AP+i. 

DEL»DEL*X/AP 


SUM-SUM+DEL 

IF(ABS(DEL).LT.ABS(SUM)*EPS)GO  TO  1 
11  CONTINUE 

PAUSE  *k  too  large,  ITMAX  too  small * 

1  GAMSER«SUM*EXP(-X+A*LOG(X)-GLN) 

RETURN 
END 
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